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ABSTRACT 


This thesis presents a new algorithm for assigning closed-loop 
eigenvalues of multivariable linear systems using output feedback con- 
trol. Subject to certain mild restrictions, the number of eigenvalues 
which can be assigned to arbitrary, distinct locations is min [m+tr-l,n] 
where m, r, and n are the dimensions of the output, control and 
state vectors, respectively. The algorithm provides an extension of 
previous results but allows a significantly large number of eigen- 
values to be assigned. 

An extensive literature survey of modal control and eigenvalue 
assignment methods based on proportional feedback control is included. 
In this survey the existing methods are classified, described, inter- 
preted and critically evaluated. Furthermore, the design options and 
design parameters in each method are identified and their usefulness 
discussed. 

Digital simulation studies involving the application of repre- 
sentative modal control methods and the new eigenvalue assignment 
method to the control of a double effect pilot plant evaporator model 
are described. In these studies the effects of the various design 
options on the response characteristics are demonstrated. 

Finally, concluding comments on the proper use of the design 
options are presented and fundamental differences between modal 


control and eigenvalue assignment are discussed. 
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CHAPTER ONE 


INTRODUCTION 


Significant theoretical advances have been made in recent 
years in the control of linear multivariable systems. However, few of 
the methods developed have been of much practical value in the unique 
control problems facing the process control engineer. First of all, 
the physichochemical characteristics of the processes are not well 
understood and thus their modeling is extremely difficult. Even if 
a reasonably accurate state space representation of the system can be 
realized and thus application of modern control theory made possible, 
the resulting model often contains a large number of interacting 
State variables of which only a limited number can be accurately 
measured. But, for the application of most of the modern control 
methods, measurement of all of the state variables is a prerequisite. 
In the relatively few multivariable control methods where Bay sub— 
set of the state variables has to be measured, the control configuration 
to be used is not specified, i.e., the design method does not provide 
any guidance concerning which variables to measure and which ones to 
manipulate. These are not the only difficulties facing the process 
control engineer, but they are the more pronounced ones. 

Modal control and eigenvalue assignment methods have been 
greatly welcomed by process control engineers since they promise to 
cope with the above problems in those cases where the state space model 
of the process is available. 


The major objective of modal control and eigenvalue assignment 
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methods is to increase the degree of stability and the speed of: response 
of the process. Modal control also provides useful guidance concerning 
the configuration of the control system which will achieve these 
objectives. 

Recently, objections have been raised against modal control 
by researchers who have developed the multivariable counterparts of 
the classical frequency domain control theory [1]. It has been claimed 
that modal control and eigenvalue assignment methods cannot deal with 
certain other control problems arising in process control. This is 
certainly true. In fact, no multivariable control method is capable 
of Boh icuane all of the desired objectives in a process control 
design problem. 

It is the author's opinion that modal control theory and 
eigenvalue assignment methods possess certain design options which can 
be purposefully used to achieve control objectives in addition to those 
of increasing the degree of stability and the speed of response of the 
controlled processes. Further developments in this area will certainly 
add to the versatility of modal control and eigenvalue assignment 


methods. 


Il “Objectives sof thesstudy 


This study is primarily concerned with modal control and 
eigenvalue assignment methods which employ proportional feedback only. 
Special emphasis is given to those methods which are applicable to 
systems with a restricted set of measurements and controls. 

The large number of publications available in this area has 
made the correlation of the available results rather difficult. In 


fact, even the distinction between modal control and eigenvalue assign- 
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ment has been largely overlooked by many. It has been one of the 
intentions of this study to classify, describe, interpret and 
critically evaluate the available methods in some detail, and thus 
prepare a basis for future work in this area. 

By combining some ideas inherent in modal control and eigen- 
value assignment ,the author has been able to devise an algorithm, which, 
subject to certain mild restrictions, significantly increases the number 
of eigenvalues assignable to arbitrary, distinct locations. It is hoped 
that this extension to the available methods will prove helpful, espe- 
cially, in its application to systems of larger dimensionality. It is 
also hoped that the basic idea employed in this algorithm will be use- 
ful in other applications. 

Workers in different fields of modern control theory have 
criticized modal control and eigenvalue assignment methods for provid- 
ing, in general, a nonunique controller matrix for the same set of 
closed-loop eigenvalues. It is one of the objectives of the simulation 
studies in this thesis to demonstrate that this lack of uniqueness may 
in fact be very useful in fulfilling other design objectives in addition 
to eigenvalue assignment. 

Another aim of this thesis is the application of some of the 
modal control and eigenvalue assignement techniques to the pilot plant 
evaporator system in the Department of Chemical Engineering of the 
University of Alberta to evaluate these techniques in this particular 
application. Although extensive simulation studies could be performed, 
the experimental evaluation of the results had to be postponed because 


of the equipment difficulties. 
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1.2 Structure of the Thesis 

The thesis has been organized in the same sequence that was 
used in the actual investigation. 

In Chapter Two, the existing methods for modal control and 
eigenvalue assignment have been surveyed, compared and critically 
evaluated. In Chapter Three, a new eigenvalue assignment algorithm 
has been presented. In Chapter Four, the results of a simulation case 
study involving the methods of Chapter Two and Chapter Three are given. 
Finally, in Chapter Five, the overall conclusions of the thesis have 


been summarized. 


1.3 Preliminary Comments 

In this thesis no attempt has been made to dwell on the 
theoretical concepts of control theory which have been extensively used 
in modal control and eigenvalue assignment methods. These can be 
found in standard texts on modern control theory, such as [2]. 

The methods described in this thesis have been based on 
the state space representation of linear, lumped-parameter, time 
invariant continuous-time systems. Their application in the simula- 
tion studies of Chapter Four has involved though discrete-time models 
derived from their continuous-time counterparts. The justification of 


this approach can be found in [3, 4]. 
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CHAPTER TWO 


MODAL CONTROL AND EIGENVALUE ASSIGNMENT: 


LITERATURE SURVEY AND THEORY 


2... jneroduction 

Modal control has been suggested by Rosenbrock [1] as a design 
technique to cope with some of the special control problems involved in 
process control. It aims at improving the stability and speed of re- 
sponse of interacting processes. Modal control is not used to eliminate 
the interactions involved in a multivariable system, instead it aims at 
controlling directly the modes of the system, and provides information 
about the control configuration that should be used for the achievement 
of its aims. 

Eigenvalue assignment techniques have been suggested as an 
approach to fulfill the same objectives. Both methods are capable of 
achieving the same objectives in cases where the designer has access to 
the values of all the state variables of the system. In all other cases 
the two methods are not expected to be equally helpful to the designer. 
Tn kate A combination of the two approaches has been successfully rea- 
lized in the control literature. 

It is the intention of this chapter to classify, describe, 
interpret, and critically evaluate these methods in some detail. The 
author hopes to identify the available design freedoms and design 
parameters in each method. The criticisms and appraisals are subjective 
and the pertinent references have been cited to help the reader to form 


his own opinions. 
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2.2 The Modes of a Multivariable System 
Consider the following state space representation of a linear 


time-invariant, lumped-parameter multivariable dynamic system, 


BCC) =e x(t) BCE Dd Cc) (2.1) 
(tie mG 2 Ct) (2. 02) 
where x(0) dis then x-1 “dimensional 6tate vector, 


uCt) "is the yr x 1° dimensional control vector, 

y(t) is the mx 1 dimensional output vector, 
and d(t) is the p x 1 dimensional disturbance vector. 

The state variables can be viewed as components of the state 
vector defined relative to the Euclidean basis of the state space. In 
general, this particular basis results in interacting state variables, 
i.e., the response of a state variable to disturbance inputs and non- 
zero initial conditions depends on the response of the other state 
variables. 


Consideration of the solution to (2.1) namely, 


A(t-T) A(t-T) 


t E 
aot = eftx(0) + f e B u(t)dt + f e= Dd(t) dz (2.3) 
O fe) 


reveals another well-known fact: Not only does the excursion of a 
controlled state variable from its steady state value introduce devi- 
ations in other controlled state variables, but also the control action 
taken to encounter this first excursion directly affects the other con- 
trolled state variables, as well. 

One possible approach to these problems, modal control, has 
been suggested by Rosenbrock [1]. It is based on the idea of finding 


another basis for the state space with the favorable property that the 
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components of the state vector with respect to this new basis are non- 
interacting. In other words, the new problem consists of the determi- 
nation of n linear combinations of the original state variables which 
are decoupled from each other. 

This decoupling basis is provided by the n linearly inde- 
pendent right eigenvectors, Wy oWo sree oWs of Ars 2), some of which 
may be generalized right eigenvectors for the case of repeated eigen- 


values [3]. The new noninteracting variables, z 72.5. ate 


y2299° 8092, 


called the canonical or modal state variables and are given by the 


linear transformation 


x= WZ rs) 
(2.5) 


or 2 = Ww 


Recognition of the fact that the right eigenvectors and left eigen- 
vectors of A, namely, Vy oVooer+ 9M» if normalized, form a biorthonormal 


Set alZdant. €s4 
<vyows> = 6), (4,5 = 1,2,-6-5) (2.6) 


allows (2.5) to be written as 


z=Vx. (J27) 


Application of this linear transformation to (2.1) and (2.2) 


yields: 


B u(t) + VD d(t) (2.8) 
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or, if A possesses distinct eigenvalues, 


mbt) wink aGhed Hut), Ed(t) (2.10) 

Varese) (2s) 

where A=V ay (2.12) 
H = ia B (2,13) 

B-V Dp (2.14) 

and F=CW (ests) 


A is a diagonal matrix, whose diagonal elements consist of the eigen- 
values of the system matrix, A. If A _ possesses repeated eigenvalues, 
then the transformation of (2.12) will either result in a diagonal or 
Jordan canonical matrix. The latter case, which involves analytical 
difficulties, may be considered as a limiting case as suggested by 
Rosenbrock [1]. 

The important outcome of the modal transformation can be 


depicted by rewriting (2.8) as: 


CAG oy ESTA ey ie Oo) eh Ca) (27.10) 
where s(t) = H u(t) (2517) 
and C8 Geel Ge) ME (2.418) 


Thus, the transformation of the control vector, u(t), the disturbance 
vector, d(t), and the output vector, y(t), into the modal inputs, 
s(t) and t(t), and the modal state vector, z(t), has created n 


decoupled first order systems, 
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from a complex system with interacting inputs and outputs [4-6]. 
The application of the transformations given in (2.4)-(2.7) 
and (2.12) to (2.3) leads to the representation of the system response 


in terms of its modal components [4]: 


n T t T -rX,T 
2 cet hy il hce, ora) ra a | es eae ar ora oath bg 
i=1 fe) 7 | 
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t 
sail v, D a(t)e : dtje 7 Wala (2.20) 
(0) 


Thus, the system response can be represented as a linear combination of 


het 
the terms e : Wis which will be called the modes of the system, i.e., 
n At n 
AGE) = Bye 6.(t) ec wo = ) 28 Cr) m.(t) C1) 
, A i : 3 al 
i=1 i=1 
ft th 
where m, (t) =e Wy is the i mode, and E(t) is defined to be 


the term in brackets in (2.20). 

This definition of the mode is not unique. Many authors have 
defined the modes of the system to be the eigenvalues or the eigenvectors 
of the system matrix, A. The definition used here is due to Zadeh and 
Desoer [2], and has been found to simplify some of the following dis- 
cussions. 

When only step changes in both the disturbance and control 


inputs are considered (2.20) can be simplified to [4]: 


T 
B Ve. Dd ASE 
x(t) =.) vy xe > - C= +(e © ) Iw, C2. 22) 


i=] eal. i 


which is possibly the most commonly encountered case in practice. 
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Furthermore, if some of the modes are complex, these modes 
contribute to the response as though they were real modes equal to the 
sum and the difference of the two complex conjugate modes. They also 
contribute sinusoidal oscillations at the frequency determined by the 
imaginary part of the complex eigenvalue [4]. 

The magnitude of the contribution made by each mode to the 
State vector is called the modal activation [4,6], and it consists of: 

i)jitecactivationm ‘of the cue mode due to nonzero initial 


states, yy x(0):s 
th 


ii)e activation of «the «i mode due to disturbance inputs, 
t » -rT 
{ v, D d(t)e sane 
) 
iii) activation of the peu mode due to control action 
t T aA, T 
2 v, B u(te dt. 
re) 


Furthermore, comparison of (2.4) and (2.21) suggests that 


ios 


xc) = z(t) Wade (2523) 


The above considerations lead to some important conclusions 
for systems with distinct eigenvalues: 
1) The modes respond to initial disturbances and process 
inputs independently of each other, i.e., to determine the modal 
; nth. ' 
activation of the i mode one makes use of the corresponding left 


eigenvector v, only. 


i 
2) Each mode is associated with a unique time constant, which 
is the negative reciprocal of the real part of the eigenvalue; for a 


system to be asymptotically stable all of its eigenvalues must have 


negative real parts. 
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3) The free response of each mode for a stable system with 
real and distinct eigenvalues follows a simple exponential decay, but 
the free response of the state variables may show a more complicated 
behavior. This is due to the fact that state variables may have 


positive components in one mode, and negative components in another. 


2.3 Ideal Modal Control 

Since the state vector is a linear combination of the right 
eigenvectors (cf., (2.23)), one way of driving it towards the origin 
of the state space (the system steady state) is by driving the time- 
varying coefficients in this linear combination, namely the modal 
State variables towards zero, i.e., by counteracting the modal acti- 
vations. This is the basic idea behind modal control. 

Thus modal control in its ideal form has to perform the 
following three tasks [6]: 

1) Determination of the value of the modal state, or equi- 
valently, the modal activations, from the system outputs (analysis of 
the outputs). 

2) Application of a linear feedback control law to the 
modal state vector which results in a modal control vector, s. This 
proportional linkage between activation and manipulation is chosen so 
that the control action is zero when the modal activation is zero. 

Thus modal control as treated here belongs to the class of proportional 
controllers. 

3)an Sy-thesiis fia control vector,:’ a,\efrom the modal control 
vector, s, such that each element of u acts on one and only one 


mode (synthesis of the controls). 
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Two comments are worth mentioning at this point: 

i) Since ideal modal control is a strictly proportional-type 
control it will invariably cause steady-state offsets in the controlled 
variables, x. 

ii) Driving the modal state variables towards the origin at 
all times does not necessarily imply that the actual state variables 
will be driven towards their steady-state values at all times. Although 
the state variables will finally reach their steady state values, there 
may be significant excursions from the origin during this process (the 
“peaking effect"). 

The action of the ideal modal controller on 2 of the system 
modes can be described by the block diagram in Figure 2.1. 

The B. andiiC .matrices in Figure 2.1 are party of the process 
and thus fixed, while the controller matrices P, K, and R may be 
chosen in any desired way, except that K must be a diagonal matrix. 

Suppose that there are n measurements (i.e., m= n) and 
n controls (i.e., r =n), and that the process is both controllable 
and observable. It is under these circumstances that all modes are 
independently available for observation and manipulation. The equi- 
valent condition which makes this possible is that both the B and 
C matrices be of rank n. 
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if th 
where P, is the i row of RP, x is the 158 column of  R;) pind 
& is any positive integer less than or equal to n [4]. Thus, if 


all the modes of the system are to be controlled the following relations 


must hold: 


Wo 

"a 
it 

i 


(2320) 


and 


les 
lbs) 

ti 
= 


(2326) 


However, if only & modes are to be controlled, then & columns of 


the 


Wg 


matrix are chosen such that the first 2 columns of BR are 
the & right eigenvectors corresponding to the 2 modes, and the & 
rows)0f §P “are chosen such tl¥at*the first = { rows of E¢ are the 
corresponding left eigenvectors. 

Under the assumption of the availability of n independent 
measurements and controls, (2.23) and (2.24) represent n equations 
in n unknowns, namely the elements of the p. and xr vectors. 
But, when the number of measurements and controls are less than n, 
it is in general impossible to satisfy (2.23) and (2.24) exactly, 
because they now represent n equations in m and r_ unknowns, 
respectively. Thus the modal activation of a particular mode cannot 
be detected, and the state vector cannot be driven in the opposite 
direction of that particular modal activation. 


The feedback control law used in Figure 2.1 is of the general 


form: 


u=-Gy. (CRP) 


Considering (2.25) and (2.26) together with Figure 2.1, (2.27) can be 


written explicitly as [6]: 


(25.8) 
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us Bal EY O -7): (2.28) 


Combining this last equation with (2.1) and (2.2), the closed-loop 


system representation in the state space is obtained: 
k(t) = (A-WKV)x(t)+Da(t) . (2.29) 
In the case where A has distinct eigenvalues, A can be written as, 
Sway (2.30) 
and (2.29) can be written as: 


k(t) = WA - K) Vv’ x(t) +Dad(t) . (2.31) 


This last equation clearly shows that the modal controller has effec= 


tively created a new system with a new set of eigenvalues, A 


das 
given by, 
dai = dj - k. (pee dal ee Tn) G2552)) 
but with the same right and left eigenvector matrices, W and a 
ine (2237) aa represents the Fee closed-loop eigenvalue and ki» 
the is diagonal element of the diagonal K matrix. 


In other words, for k. > 0, the time constants associated 
with the modes have been reduced, but the directions of the modes have 
remained unchanged. This is the net effect of the ideal modal control- 
lers on the system eigenproperties. 

The effect of the ideal modal controllers on the system 
response can be demonstrated by considering a step change of inputs 
at time zero. For the case where & of the n _ modes are controlled 


the response of the closed-loop system can be obtained from (2.22) [4]: 
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Obviously, the factors (pt in the first summation reduce the ae 


1 EL 


activation to a great extent for large values of ke Thus, for inputs 
driving the process away from steady state, the ideal modal controller 
wild decrease the magnitude of the undesirable contributions of the 2% 
controlled modes to the system response. 

Another interesting property of the modal controllers can 


be depicted by rewriting (2.31) as [6]: 


u 
x(t) =Ax- ) kalsve ,wkehwiattdDed (2334) 
i=1 
x 
or x(t) =Ax- ) kuthu xiteDed (2435) 
= et 5 ae er ae 
i=1 

where Fi = Wee par denotes the dyadic product of the right and left 
eigenvectors corresponding to the ‘bee mode. Thus, ideal modal con- 


trollers represent the sum of a set of dyadic single mode controllers, 


namely, ue ) k EB x. Associated with each mode to be controlled 


there is a meaSurement vector, a manipulation vector, w and 


ae? oe? 
a gain Kis which directly reflects the amount by which the correspond- 
ing open-loop eigenvalue has been changed. 

The modal control, through the use of modal analysis and 


synthesis, can be employed to shift each open-loop system eigenvalue 


indpendently to a prespecified closed-loop location. Eigenvalue assign- 
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ment has a similar aim: replacement of the open-loop eigenvalue set 
with a more desirable set for the closed-loop system. In those instances 
where the designer has access to all the values of the system state 
variables, modal analysis has been employed in eigenvalue assignment 
to gain control over the fate of each individual eigenvalue. 

Modal control changes the system eigenvalues by controlling 
its modes while eigenvalue assignment changes the system eigenvalues 


without directly controlling its modes. 


2.4 Use of Ideal Modal Analysis in Eigenvalue Assignment 

The control method to be described in this section is mainly 
due to Simon and Mitter [7]. It performs two of the tasks of ideal 
modal control. It determines the modal state variables and uses feed- 
back to assign the system eigenvalues without directing the controls 
along the right eigenvectors. This has several implications: 

1) Some or all of the system eigenvalues can be positioned 
without changing the remaining ones. The necessary and sufficient 
condition required is the complete controllability of the modes whose 
eigenvalues are to positioned. 

2) The modes will be interacting and at least some of the 
eigenvectors of the system will change direction arbitrarily. 

3) The control energy required for a given eigenvalue shift 
is greater than that required by an ideal modal controller [8]. 

4) Measurement of all of the state variables is necessary, 
but the number of controls, r, can be less than the number of states. 

The theoretical basis of modal analysis and eigenvalue assign- 


ment has been presented by Simon and Mitter [7], Wonham [9], and Mitter 


qiotem af nokisse sidt at lindbsobb od ox/boritem Lozsneo sdf 

fssbi to s¥ap3 of} to ows dapettem 21. . [1], 76d9%H bas aomi2 02) ab 
-b692 seen bas sefdelrey state Ishom sii senimyetsh 32 .foxsn0o Iabom 
aloxines of3 gnttsstib iupdtiw avulevasgie meteye oda aglees of doed 
:enoitestiqmt isreyee esi etd .#rossovaegts tigtx odt gnole 
benettieog od aso esulawnegis mezeye ait to iis xo amoe (i a) 
juoits b¥tue. Gls Yraese0en sft .eene guiniemey old gaigasds aiodjiw 
seodw eébom 4it to winner etofqmos on? al bartepex motsthnaos 


benotztaog 0} 918 esulavasgis. 


ads Yo aos tunel to Bae gultonrsim od 1thw edbom adt (© 


18 


and Foulkes [10]. Two extremely useful design tools in the eigenvalue 
assignment problem have been introduced by Simon and Mitter [7]; the 
dyadic controller and the recursive eaten technique. Since these tools 
have been extensively used in subsequent studies, Simon and Mitter's 
work [7] deserves careful consideration. In the following discussion 
it will be assumed that the open-loop system has distinct eigenvalues 
or has been converted into one with distinct eigenvalues as described 
aha dy) bib BA i 
Three basic design steps are involved in Simon and Mitter's 
work [7]: 
i) Conversion of the system equations into their modal 
canonical form via a similarity transformation. 
ii) Construction of the constant controller matrix providing 
for the desired eigenvalue movements in the modal apace. 
iii) Conversion of this controller into the form u = G Xx. 
The system representation obtained through the application 
of the first step will be that given by (2.10). Next, define the 
control law to shift & of the system eigenvalues (2 <n) by feeding 


back the corresponding %& modal state variables 


ti = 


Hq 2 


(2536) 


or equivalently 


~ ~ 


u= 8) 2, + g Zo = es By 2 (23.354) 


where Be ahye ties are r-dimensional vectors representing the first 
& columns of GQ. The remaining (n-%) columns consist of O-vectors. 


The resulting closed-loop system in the modal domain is given 


by 
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or equivalently [7]: 
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Equation (2.39) indicates that the last (n-%) eigenvalues 
of the system are unaffected by the selected control scheme. On the 
other hand, the first & eigenvalues can be assigned by choosing 
By Boor By such that the upper left partition of the system matrix 
in (2.39), denoted by A » has the desired eigenvalues. This can be 


achieved by the following 3 steps: 


Nt 


1) Find the characteristic equation of 


dett [ai 7— A] = a +a oe +O ota. 28a =O 2 (2.4)) 
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Note that (2.39) and (2.41) imply that, in general, ay Clie eel 2 Fe eee) 
are nonlinear functions of the components of ee Clem Lele sae a Ut, 
in the important special case where r=1, they are linear functions 
f Prat | 
of the 8, 8. 
2) Obtain the characteristic equation corresponding to the 


}: 


desired eigenvalue set, {Aa hao? ae ve d ap 


Wea ie eee (2.42) 


(S-A 44) (S~Ago) «+ (SHA fi f-1 8 


dd 


3) Compare the two characteristic equations of Steps 1 and 
2. For the closed-loop system to have the desired eigenvalues, the 
following condition must hold: 


EOE spend Un EI V Saree tes Sy ee (2.43) 


Equation (2.43) consists of a set of & nonlinear equations 


in the (+r) unknown components of 


WGr 


Dyadic Approach: 


The above discussion reveals the difficult computational 
problem involved in the method and suggests a possible solution to it, 
namely conversion of the multi-input system to an equivalent single- 
input system in order to convert a nonlinear algebraic problem to a 
linear one. This idea, which forms the basis of the dyadic approach, 
is realized by restricting the controller to have linearly dependent 


columns. Thus the suggested new control law has the form [7]: 
uu = [kj 8, Kk,» aha. k glz) (2.44) 


or simply 


x (2.45) 


($8.8) St) gd. es ‘a - eke: 


tad sain dete algae a 
of etusavingtin tor tesh witi ovadd 60 setige qool-Aaints wi ~ roe oe 


& ch i ges :blod teum cotitbood gatwolfot 


(€a.$) » (hg oo ng Sd dom, 2) +4 * ls | ih ie e 7 


¢ 


eaotisups teentiaon 4 to 392 s to extetanos (&8,S) noksaupa | 


cs 1 Yo eanonognas awonslau aes ot ak 
7 


-slonoxggh 9tbeya 
fsnottssuquos tluolttib eda elasvex moteausetb, svods "edt an 


eit o3 mottulos sidteaog s aseoggue ~~ bodsem et at bevioval ; 
> * they »e : = 7 
-signts 2nafeviups on 09 moseve tugni=t3Lun ae 30 notexsvno vloasa iu 


“+ 02 moldong olasdngin saan » ayovaeo 03 zobxe of mate 


ij s © oh %Ave - “ke 


.onorage atbeeb ont 20 wtnad as . oki yeobt abet pet 
he 2 a ® =ey : 


weak ae Lf 
; f mii 


Biers 


: 


21 


if it is desired to assign & of the system eigenvalues corresponding 

; ft ‘ 
to the & modal state variables, z.. Yo contains the & left eigen- 
vectors corresponding to the eigenvalues to be shifted, as its 2 rows. 


The equivalent single input system to that expressed by (2.1) 


is thus given as: 


0) Rex Ct) buGe) |r oDva (1) (2.46) 

where 
be Bg (2.47) 
and u(t) =k Vx . (2.48) 


The choice of g is arbitrary as long as the equivalent 


system in (2.46) remains controllable, i.e., 
eer WRN eR ebay iat (2.49) 
iS 9-9 Og ° “ 


The choice of Ki, on the other hand, is dictated by the desired set 

of eigenvalues. This problem has been considered by Simon and Mitter 
[7], Gould, Murphy and Berkman [12], and Retallack and MacFarlane [13]. 
The solution, which seems to be the most straightforward and computa- 
tionally easy one [13], is based on a result of Hsu and Chen [14] 
relating the ratio of the closed- and open-loop characteristic equations 
to the determinant of the system return-difference matrix. 


The elements of x are given by [13]: 
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Design Freedom: 


The arbitrariness in the choice of g can be utilized as 
a degree of freedom for one of the following design objectives: 

1) eliminating the need for certain controls, 

2) eliminating the need for certain measurements, 

3) maximizing the "measure of controllability" of the shifted 
eigenvalues thus minimizing the absolute value of the controller gains 
to be used, 

4) affecting the transient response of the system. 

This first objective can be realized by simply assigning 
the value of zero to the elements of g corresponding to the control 
elements to be neglected providing that the controllability condition 
given by (2.49) is not violated. 

Simon and Mitter [7] have suggested that (r-l) of the n 
state measurements can be deleted without affecting the attainability 
of n eigenvalue locations. This can be achieved by choosing the 
gain elements in x! such that the row vector a v,) defined in 
(2.48) has (r-1) zero elements. These zero elements will correspond 
to those state variables, which now do not have to be measured. Then 
g can be determined to assign the desired eigenvalue set. 

The third objective is of great practical importance, since 
high controller gains may lead to the unsatisfactory control of a real pro- 


cess due to noise effects, model inaccuracies, nonlinearities, and satur- 
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ation of manipulated variables. In the event of shifting one system 


eigenvalue, say A by using one control element only, a measure of 


ik: 
the degree of controllability of Ay by the scalar control u_ con- 
sidered is given by [7]: 
even Be 
V4 S| fog eae eae (2553) 
Evie Sqn 


where vy is the left eigenvector corresponding to ST and b is 
the control coefficient vector corresponding to the control element u. 


V5 has been generalized by Simon and Mitter [7] to the case of dyadic 


control using all the control elements by: 
vy ——————— m2 54) 


where [| gl | - max ls, | by definition. 


Siyaiwht 
The above measure of controllability is maximized through 


the choice of g as [7]: 
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The concepts of the degree of controllability and the degree 
of observability have been considered for multi-input multi-output 
systems in [15-18]. But the problem of adjusting the variable system 
parameters, e.g., g, such that these degrees are maximized seems to 
offer challenging mathematical difficulties [17, 18]. 

The fourth point is also of great stati oa importance. 
Several authors [19-22] have already noted that for multi-input systems 


the controller matrix assigning a given set of eigenvalues is non-unique. 
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In particular each choice of g will result in a different controller 
matrix. Thus, in general the eigenvectors of the closed-loop system 
will depend on the choice of g which implies that for the same set 

of closed-loop eigenvalues, different transient reponses will be 
obtained. The analytical relation between the direction of g and 

the closed-loop eigenvector directions, and thus the shape of the tran- 
sient response, is not known. 

The dyadic approach above described simplifies the computa- 
tional aspect of eigenvalue assignment, but this advantage is gained 
at the cost of losing some design freedom in specifying G P1Sktx2 14 
The recursive design technique of Simon and Mitter [7] offers a solution 
to this problem and also reduces the computational effort if only a 
subset of eigenvalues are to be assigned. 

The recursive dyadic approach of Simon and Mitter [7] ree 
of a series of steps where at each step a small subset of eigenvalues 
are shifted by the single-step dyadic approach and the closed-loop 
system matrix resulting from the last step is considered to be the new 
open-loop system matrix. Assignment of & eigenvalues by this approach, 


one at a time, will lead to the following control law: 


Me hy ie ee (2.56) 
where u, = k vi x (ho2t 52 oF. 52) (2757) 
SB i a= se a | 
r - i . 
ae ey) ee ee (Tie dl pied csi (2.58) 
<B Vi» £.7 


The determination of g. to satisfy some design objectives in addition 
to eigenvalue assignment is an unresolved problem in control systems 


theory. 
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Simon and Mitter's [7] work has been very useful from a 
theoretical point of view, and has provided the basis for more practical 
methods, which attempt to assign only a small subset of system eigen- 
values and require only a small number of measurements. These methods 


are discussed in Sections 2.5, 2.7-2.10, and in Chapter Three. 


2.5 When the Number of Measurements and the Number of Controls are 
Less than the Number of States 

The application of ideal modal control requires the avail- 
ability of n independent measurements and n independent controls, 
but in industrial situations this requirement usually cannot be ful- 
filled due to either physical or economical considerations. Thus, the 
problem of conditioning the system modes in the absence of n indepen- 
dent measurements and controls is an important one. 

Measurement of a small number of state variables leads to the 
projection of the state vector onto a subspace of the state space, 
called the measurement space. To determine the modal states, one must 
be able to measure the components of the state vector along the right 
eigenvectors (cf., (2.4), (2.5)). This is done by determining the pro- 
jection of the state vector along the corresponding right eigenvector. 
But the deficiency in dimension of the measurement space prohibits this 
and leads one to determine the projections of the right eigenvectors 
onto the measurement space instead. These projections are then used 
in place of the right eigenvectors. The resulting set of vectors, 

(wis Wh s mA we) is linearly dependent. Consequently, the projection 
of the left eigenvector of the 1B controlled mode onto the measure- 


ment space, namely v* will not be in general orthogonal to the set, 
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Ad wi. Thus, measuring the activation of a 
controlled mode will, in general, result in a quantity that contains 
the activation of all the modes. 

The activation of a mode can be estimated reasonably well 
by the activation of the mapping of the mode in the measurement space 
if the components of the mode that are orthogonal to the measurement 
Space are small in magnitude compared to the mode itself. If any mode 
is orthogonal to the measurement space, then that mode is completely 
unobservable. Therefore, the measurement space must be chosen such 
that the mode whose activation is to be measured has small orthogonal 
components. 

In ideal modal control if only one mode is disturbed, then 
the resultant control action is directed along the appropriate right 
eigenvector, so that the component added to the state velocity vector, 
aGens directly counters the disturbance. Thus, the control action 
does not influence the other modes. This cannot be achieved if 
thetrankhibh=wri< myesince then the manipulation space spanned by 
the columns of B cannot, in general, contain the appropriate right 
eigenvector. In fact, generally all of the right sastinastis will be 
outside of the manipulation space. It is also worth noting that, if 
a mode is orthogonal to the manipulation space then that mode is 
uncontrollable. 

In summary, the small number of measured variables means 
that estimates of modal activations cannot be the true activations, 
and the small number of manipulated variables means that the control 
action cannot exactly counteract the activations of the controlled 


modes. As a result the control action produces disturbances in modes 
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other than the controlled mode [4, 6]. 

Two basically different but practically quite similar approaches 
have been adopted to cope with these difficulties, approximate modal 
control and eigenvalue assignment techniques. 

The major reasons for the inability to achieve ideal modal 
control in the absence of n independent measurements and controls 
are the difficulties involved in estimating a specific modal activation 
and in influencing a specific mode without causing any modal interactions. 
Consequently, there is considerable motivation for trying to alleviate 
these difficulties in order to achieve some of the favorable results 
obtainable from ideal modal control. This approach will be classified 
here as "approximate modal control". 

The major effect of modal control, which accounts for its 
success, is the shift of the open-loop eigenvalues to more favorable 
locations in the closed-loop system. Thus, one should try to change 
the system eigenvalues disregarding any interaction between the modal 
states which will result from this practice. This approach will be 
classified here as eigenvalue assignment. Eigenvalue assignment tech- 
niques simply insure that some set of desired eigenvalues form a sub- 
set of the closed-loop system eigenvalues. Thus, they do not provide 
any pairing between open-loop and closed-loop eigenvalues, and cannot 
control the movement of the unassigned set of system eigenvalues. 

Nor do they attempt to preserve the original eigenvectors. 

Eigenvalue assignment methods using approximate modal 
analysis try to pair open-loop and closed-loop eigenvalues and to 
maintain the unassigned set of eigenvalues at their open-loop locations. 


Several successful simulation studies employing approximate 
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modal control have been reported. This is mainly due to the fact that 
modal control provides the designer with criteria for deciding which 
modes to control and which measurements and manipulated variables to 
use for this purpose. The success of the resulting modal controller 
largely depends how well one can meet these criteria. Eigenvalue 
assignment, on the other hand, does not provide the designer with 
direct information related to the control configuration that should 
be used. But, it is possible to combine modal analysis with eigen- 
value assignment, and this in fact leads to what has been classified 
as "approximate modal analysis" in the preceding paragraph. Approxi- 
mate modal analysis as developed by E.J. Davison and his co-workers 
[23-25] is, in the author's opinion, one of the more practical and 
powerful techniques described in this thesis. There are, of course, 
special applications in which some of the other techniques will 


prove to be more useful. 


2.6 Criteria for the Selection of Dominant Modes, Measurements, 


and Controls 

Physical systems with a few well-separated dominant modes 
have been reported to perform better under modal control and eigenvalue 
assignment schemes using modal analysis than under conventional 
control schemes [1, 4, 6, 23-28]. Thus, the designer should first 
analyze the system modes to determine the dominant ones, namely, the 
ones which affect the system response the most. This is extremely 
important in view of the fact that the control of every additional 
mode results in much greater difficulties in the design stage. 


Consideration of (2.20) indicates that there are three 
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attributes associated with a mode which must be analyzed: 

1) The right eigenvector corresponding to the mode. Com-~ 
parison of the ve element of a normalized right eigenvector to the 
hb element of the other right eigenvectors indicates the relative 
influence of that mode on the ie State variable [29]. 

2) The eigenvalue of the mode. This indicates the time 
constant associated with the mode. 

3) The modal activations. These indicate the influence of 
a mode on a state variable for specific disturbances and initial con- 
ditions. 

The first two attributes are the more important ones to con- 
sider. It is desirable to control those modes which have the largest 
influence on the important process variables and among these, the 
modes with the larger time constants. In case some of the modes 
are unstable, they must definitely be included among those to be 
controlled. 

Some typical processes have a wide spectrum of eigenvalues, 
thus the choice of the modes to be controlled might seem easy. But 
it is quite possible that the occurence of a disturbance will create 
dominant modes out of less important ones by highly activating them. 

The next question which needs some consideration involves 
the choice of the "best'' measurements and manipulated variables 
from the often restricted set of possible ones in order to control the 
dominant modes. 

The "best'"' measurements, in the modal sense, are those which 
are most sensitive to the controlled modes and least sensitive to the 


uncontrolled modes. The "best" controls are those which have the most 
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effect on the controlled modes while having the least effect on the 
uncontrolled ones. 

The above definitions are applicable only if the measurement 
and control sections of the modal controller are noninteractive. This 
being true for the ideal modal controllers only, the following criteria 
for the choice of the measurement and control spaces are only approxi- 
mately true; 

1) The measurement space should be chosen such that com- 
ponents of the controlled modes which are orthogonal to it are of 
minimal magnitude and the projections of the uncontrolled modes onto 
this space have minimal length. 

2) The control space should be chosen such that the pro- 
jections of its basis vectors on the controlled modes are of maximal 
length while being of minimal length on the uncontrolled modes. 

These arguments are another manifestation of the desirability 
of forming measurement and control sections, which in the absence of 
n independent measurements and controls, act like the left and right 


eigenvector sets of the original system. 


2.7 Approximate Modal Control 

Approximate modal control aims at changing the dominant 
eigenvalues of a system, which has fewer than n measurements and 
n controls, without changing any of its other eigenproperties. 
Basically, the problem consists of designing measurement, P C and 
control, BR, sections, which in the absence of n measurements and 
n controls, play the same role as their counterparts do in an ideal 


modal controller described in Section 2.3. Two different approaches 
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have been used in attacking this problem: 

1) On the basis of the criteria given in Section 2.6, choose 
the "best" measurement and control variables, thus fixing the C and 
B matrices. Then, determine the P and R matrices in one of the 
alternative ways to be described later. 

2) Rather than first choosing the "best" measurement and 
control variables and then determining the measurement and control 
sections, perform these two tasks simultaneously. 

The first approach has been adopted by Rosenbrock [1], and 
by Takahashi and his co-workers [5, 30]. In describing these methods 
it will be assumed for simplicity that the number of measurements and 
controls are equal (m= r = &) and the open-loop eigenvalues are 
distinct, the generalization of the methods for the case where m#r 
is straightforward and can be found in the original works [1, 5, 30]. 

Rosenbrock [1] suggests choosing the P matrix in 


Figure 2.1 as 


) C2259) 


where Wo contains as its columns the 2 right eigenvectors cor- 


responding to the & dominant modes of the system. Similarly, the 


R matrix in Figure 2.1 should be determined from 


it = 
& 


R= (V) B) (2.60) 


litt 


where v, contains as its rows the & left eigenvectors belonging 


to the & dominant modes of the system. Thus, the approximate modal 


controller designed by Rosenbrock's method becomes: 
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(2561) 


where K is an & x & diagonal matrix, whose diagonal elements are 
to be obtained from (2.32). This method has been applied in simulation 
studies to the control of distillation columns [1, 4, 26, 28] and to 
the control of a boiler [27]. In all of these studies it has been 
observed that the control of a very small number of dominant modes 
gave a Significant improvement over that obtainable by corresponding 
conventional controllers. Levy's work [4], which is mainly concerned 
with the modal analysis of different distillation column models, is 
extremely useful for any further work in this area. It effectively 
demonstrates the usefulness of modal analysis in determining the 
control configuration of a system. 

Takahashi et al [30] have realized that the objectives of 
approximate modal control as defined at the beginning of this chapter 
could be realized if the (n-%) uncontrolled modes of the system were 
uncontrollable and unobservable while the remaining %&% modes are both 
controllable and observable. For the type of systems under consider- 
ation this would imply that the H and F matrices defined by (2.13) 


and (2.15) are of the following special form: 


Ho H 

a lee Gaee bean (2.62) 
ea: y 

Bm siReiiaie l= Spy, wild] de: (2.63) 


Any system possessing this characteristic has been designated by 


Takahashi et al [28] as one with "ideal measurement and control" sections. 
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The corresponding modal controller is again the one given by (2.61). 
If (2.62) and (2.63) are fulfilled, this controller performs the 
functions of an ideal modal controller for the 2% modes of the system 


provided that 


rank Hy = rank Ey = s CZa04) 


If, on the other hand, only (2.62) and (2.64), or (2.63) and (2.64) 
hold, then the values of & of the system eigenvalues can be changed 
without changing the values of the remaining eigenvalues, but nothing 
can be said about the directions of the closed-loop eigenvectors. 

The conditions for ideal control and measurement sections 
as expressed by (2.62) and (2.63) are very difficult to satisfy. The 
failure to fulfill either (2.62) or (2.63) implies that the controller 
of (2.61) cannot exactly fulfill any of the objectives of an ideal 
modal controller. This can be demonstrated by considering the modal 


domain closed-loop system matrix: 


: -1 
A \ ae Oy BA 
Re oti ae ea per ionsy re ae (2.65) 
vl az = 
- 4 HY K =n-2 Au Ho K Fo AF 
where A is an & x & diagonal matrix with diagonal entries con- 


=Ld 


taining the & desired eigenvalues of the closed-loop system. A 
is an (n-%) x (n-%) matrix containing the (n-%) open-loop eigen- 
values corresponding to the uncontrolled modes. Thus, the performance 
of the controller of equation (2.61) will depend on the designer's 
ability to find such measurements and controls which minimize the 


elements of AF and AH. 
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Loscutoff [5] has suggested specifying nonzero off-diagonal 
elements in K such that either K Fé AF oor 0H Hu K would become 
a null matrix without affecting the assignability of 2 eigenvalues. 
This approach necessitates the solution of nonlinear algebraic equations. 
Also, it cannot be employed unless the condition, mtr > ntl, holds. 
Whenever this last inequality holds all of the n_ system eigenvalues 
can be assigned. This approach does not provide any improvement over 
Simon and Mitter's [7] suggestion to reduce the number of state measure- 
ments (see Section 2.4), in cases where the state variables rather than 
the outputs are measured. 

Apart from their contribution to modal analysis, Takahashi 
et al [30] have suggested the derivation of another approximate modal 


controller, whose analyzer and synthesizer are obtained from (2.25) 


and (2.26) by using right and left pseudoinverses. Thus, 
B= P(E F) (2.66) 


were (2.67) 


ERCrer) (2.68) 


where K is an nxn _ diagonal matrix, whose diagonal elements are 
obtained from (2.32) of this chapter. The existence of the pseudo- 
inverses in equations (2.66) and (2.67) depend on the following rank 


conditions 


pankolu=crankole=e2.e C2469) 


Takahashi et al [30] have not supplied any physical or mathematical 
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conditions under which this method would perform satisfactorily. This 
is the major weakness of the method; one does not have any a priori 
knowledge about the expected performance of the controller. 

The second approach mentioned at the beginning of this 
chapter has been developed by Howarth [6], who has also applied the 
method to control one or two modes of the distillation column models 
developed in [4] and [28]. 

Howarth has defined two objective functions, which give a 
quantitative measure of the functional resemblance of the measurement 
and control sections to be synthesized to those of an ideal modal 


controller. The objective function defined for the measurement section 


is 
<a pi w.> 
¢ = max | —————— Gipsy 252 nD (2270) 
j#i <c! . w.> 
= Pie i 


B = max (tad Pe ye) C2re 7s) 


4#i 


Hod | ltd 
tn) 
Vv 


where i refers to the controlled mode, P. to the ice row of P, 
x to the as column of R, where P and R have been defined 
in Figure 2.1. Minimization of a with respect to the possible 
measurements selectable from a fixed total set and with peepee to 
the elements of B. fixes the measurement section of the approximate 
modal controller for the ey mode. Similarly, minimization of 8 


with respect to the possible controls selectable from a fixed total 


set and with respect to the elements of x, fixes the control section 
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of the approximate modal controller for the tan mode. Howarth [6] 


Suggests attacking the multiple-mode control problem as a series of 
independent single-mode control problems, which can be solved by the 
above procedure. This is based on the assumption that the various 
single mode controllers thus designed will not interact in the modal 
denaurt Apart from this the following assumptions have been made in [6]: 

i) The measurement and control sections are noninteracting. 

ii) The eigenvectors of the closed-loop system are the same 
as those of the open-loop system. 

Howarth [6] has drawn several conclusions from his simulation 
studies involving the application of his technique to distillation 
column control: 

1) For the case of single-mode control none of the eigen- 
values of the system move towards less desirable locations. Thus, 
regardless of the gain values used, the natural stability of the system 
is preserved. 

2) The controller approaches the performance of the ideal 
modal controller at low gain values. 

3) The single-mode controller moves the zeros of the system 
close to the eigenvalues of the uncontrolled modes. This implies high 
sensitivity of the closed-loop system to parameter variations. 

4) For single-mode controllers the eigenvalue of the con- 
trolled mode approaches that of the next slowest one as the value of 
the gain element is increased. 

5) To approach the performance of an ideal modal controller 
for control of a single-mode, the eigenvalue of the slowest mode must 


be well separated from the other eigenvalues, and the eigenvectors must 
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be well distributed in the state space. 

6) The performance of multiple-mode controllers sysnthesized 
by this approach is, generally, considereably less satisfactory than 
that of single-mode controllers. 

Among the approximate modal controllers the one proposed by 
Howarth [6] comes closest in its performance to the ideal modal con- 
troller, but it is rather difficult to justify the computational efforts 
required in the design process. In fact, any modal controller which 
Succeeds in reducing the time constants of the slow modes of the 
system and does not create any new slow modes may very well be worth 
considering. The identity of open- and closed-loop eigenvectors, which 
among the approximate controllers, is best achieved by Howarth's 
method, is a rather unjustified design criterion. It is possible for 
a controller to result in closed-loop eigenvectors which insure a 
better performance than the open-loop ones, e.g., better disturbance 
rejection properties, lower sensitivity to controller parameters, better 


steady state performance, etc. 


2.8 Use of Approximate Modal Analysis in Eigenvalue Assignment 


Davison and Goldberg [23] have combined modal analysis with 
the output feedback eigenvalue assignment technique of Davison [31] 
to create the output feedback counterpart of Simon and Mitter's [7] 
technique. It has been rather successfully applied in simulation 
studies by Davison and his co-workers [23-25]. 

Modal analysis involves determination of the modal state 
vector which, in general, involves the measurement of all of the state 


variables. Once some of the modal state variables are known, their 
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proper feedback will ensure the assignment. of the corresponding eigen- 
values without disturbing the remaining ones. Thus, there is consider- 
able incentive to determine the modal state variables. Davison and 
Goldberg's method [23] involves generating an approximation of the 
modal state variables via the measurement of a few state variables. 
This may be achieved by selecting those measurements which at normal 
operating conditions contribute most to the particular modal states 
whose eigenvalues are to be changed [23]. Now, feedback of m modal 
States, where m is the number of measurements, by use of the output 
feedback law devised in [31] will guarantee exact assignment of the 
m corresponding eigenvalues. Thus, if the modal state variables 
corresponding to the dominant modes can be approximated well enough, 
then the dominant eigenvalues will be changed. But, since in the 
absence of n measurements the modal states cannot be estimated 
exactly, the remaining (n-m) eigenvalues will move from their open- 
loop locations, too. 

The successful application of approximate modal analysis in 
eigenvalue assignment depends on how well the dominant eigenvalues 
are separated from the rest and how well one can estimate the modal 


states. 


2.9 Exact and Approximate Eigenvalue Assignment via the Direct Use 


of the Process Measurements 
Eigenvalue assignment techniques involve synthesis of a 
closed-loop system which possesses a specified set of eigenvalues. 
Unlike modal techniques, eigenvalue assignment. techniques do not 
attempt to attribute a physical significance to each eigenvalue and 


are not directly concerned with the change in eigenvector directions. 
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Most of the well-known eigenvalue assignment techniques have 
adopted the objective of exactly attaining a specified closed-loop 
eigenvalue set. This objective seems to be an unnecessarily ambitious 
one at this time. The behavior of the closed-loop system depends on 
both its eigenvalue locations and its eigenvector directions, as 
indicated by (2.20). But, currently there is no systematic way of 
selecting these in order to guarantee certain performance characteristics. 
Also, the problem of synthesizing closed-loop eigenvectors has neither 
been solved nor attempted in the control literature. 

A less ambitious objective employed in other eigenvalue 
assignment techniques is that of shifting the eigenvalues to some region 
in the complex plane rather than to specific locations. This region 
might be the left half of the complex plane or some region around a 
particular point. Techniques which try to achieve this objective will 
be referred to here as "approximate eigenyalue assignment techniques”. 

The process measurements are commonly all of the state vari- 
ables, some subset of them or some functions of ee Feedback control 
techniques employing measurements of all of the state variables or n 
independent functions of them (i.e., the rank of matrix Gy ain (232) 
is n, which in turn is the order of the system) are termed as state 
feedback techniques. Feedback control techniques employing some sub- 
set of the state variables are termed incomplete state feedback, and 
those employing m (m <n) independent functions of them (i.e., 
rank € = m) are termed output feedback techniques. Here, no dis- 
tinction will be made between incomplete state and output feedback 
techniques, and they shall be referred to as output feedback techniques. 


Most of the chemical processes which present control problems 
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can be represented by high order ee It is also the rule, rather 
than the exception, that some of the state variables involved in these 
models cannot be measured directly with a desired amount of accuracy. 
In fact, economical considerations often put a stringent limit on the 
number of state variables that are measured. Thus, methods based on 
the assumption of availability of all the state variable measurements 
are not often of direct practical utility. Two alternatives available 
to the designer are: 

1) Output feedback methods based on the original high order 
model. 

2) State feedback methods based on a reduced order model. 

Output feedback methods capable of handling both exact and 
approximate eigenvalue assignment have been devised. They are capable 
of fulfilling the eigenvalue assignment objective in cases where 
access to all of the state variable measurements is possible. 

Model reduction techniques which are most appropriate for 
both modal control and eigenvalue assignment’ are those based on the 
modal characteristics of the system. Representative methods in this 
category can be found in [32-34] and numerous related articles. The 
success of model reduction in a particular system depends on how well 
the dominant modes retained in the reduced model are separated from 
the neglected modes. Use of model reduction in conjunction with modal 
control and exact eigenvalue assignment will almost invariably result 
in the approximate rather than exact assignment of the system eigen- 
values. 

The number of eigenvalues which can be assigned exactly in 


a particular system crucially depends on the number of available 
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measurements and controls relative to the order of the model and on 
the number of controllable and observable modes. State feedback tech- 
niques guarantee the assignment of all those eigenvalues that belong 
to controllable modes [7, 9]. For output feedback techniques, on the 
other hand, the necessary and sufficient conditions for eigenvalue 
assignment are more stringent. Davison [31] and Davison and Chatterjee 
[35] have provided the sufficient conditions for the assignment of m 
and max (m,r) eigenvalues, respectively, where m and r_ are the 
number of nontrivial measurements and controls. The sufficient con- 
ditions for the assignment. of (mtr-1) eigenvalues have been 
developed as part of this thesis and can be found in Chapter Three. 
The more practical state feedback algorithms for eigenvalue 
assignment have been first developed by Simon and Mitter [7], Porter 
and co-workers [36-39], and by Anderson and Luenberger [40]. The 
algorithms in [7] and [36-39], which are very similar, have been 
derived from modal analysis. Their description has been attempted in 
Section 2.4. The extension of these methods to systems with complex 
eigenvalues can be found in [38-39]. Anderson and Luenberger's method 
[40] is based on the transformation of the system matrix to its phase 
variable canonical form. The linear transformation leading to the 
phase variable form (unlike the modal transformation) is nonunique 
for multi-input systems. This fact may be utilized as a design free- 
dom [41, 42]. Also, this transformtion and the associated eigenvalue 
assignment methods do not necessitate the determination of the system 
eigenvalues and eigenvectors. These are a few advantages of Anderson 
and Luenberger's [40] method, but it has also disadvantages: The 


computational effort is almost independent of the number of eigenvalues 
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to be placed, the system must be completely controllable or else 
decomposed via Kalman's structural theorem [43] before this eigenvalue 
assignment technique can be applied (this difficulty can be avoided 

by the use of the canonical form suggested in [44]), matrix inversions 
are required, and it does not possess the clarity and simplicity of 
Simon and Mitter's method [7] and its later modifications [12, 13]. 

The feedback matrix guaranteeing the assignment of a set 
of eigenvalues is not unique in multi-input systems. This design free- 
dom can be exploited to fulfill other design criteria than assignment 
of a set of eigenvalues. In Simon's methods the choice of the g 
vector in (2.45), the pairing of open-loop and closed-loop eigenvalues, 
the number of eigenvalues shifted at each step of the recursive design, 
the sequence in which the open-loop eigenvalues are shifted and the 
number of control elements used at each step of the recursive design 
technique constitute. the major design decisions. These design free- 
doms can be very important in obtaining a "better" transient and 
steady state response, a lower sensitivity to parameter variations, 
and a higher degree of integrity. But, it is unfortunate that systematic 
methods are currently not available to exploit these design freedoms, 
i.e., the designer has to tune the controller. 

The more well-known output feedback techniques devised for 
exact eigenvalue assignment are due to Davison [31], Davison and 
Chatterjeee [35], and Jameson [45]. Sridhar and Lindorff [46] have 
provided an alternative constructive proof to Davison's theorem [31], 
but the validity of their proof has been questioned by Topaloglu and 
Seborg [47]. 


The techniques of [31, 45] ensure that the open-loop set of 
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eigenvalues is replaced by a closed-loop set some of whose elements, 
the assigned eigenvalues, are prespecified while the remaining eigen- 
values cannot be arbitrarily assigned. Experience with chemical pro- 
cess systems [4, 6, 23-28] indicates that changing the value of a few 
"badly'' located dominant eigenvalues of a system can improve its per- 
formance significantly unless the rest of the eigenvalues move to 
"highly undesirable" locations. The eigenvalue assignment methods 

of [31, 35, 45] may not always fulfill these conditions since they do 
not ensure control over individual eigenvalues and cannot dictate the 
fate of the unassigned eigenvalues. Thus, it is very possible for a 
designer to end up with a closed-loop eigenvalue set which is worse 
than the open-loop one, because he has inadvertently moved the well- 
located eigenvalues to better positions and the "poorly" located eigen- 
values being left uncontrolled have moved to worse locations. 

The methods of [31, 35] and [45] can be improved by modifying 
them to be able to gain control over the individual eigenvalues, to 
increase the number of exactly assigned eigenvalues and to approximately 
assign a4 rest of the eigenvalues. The first objective may be achieved 
to some extent by the use of approximate modal analysis as described 
in Section 2.8. In cases where this approach can be applied satisfac- 
torily, the controlled eigenvalues will be shifted to desirable closed- 
loop locations without shifting the uncontrolled eigenvalues. The 
second objective has been achieved to some extent through the algorithm 
developed as part of this thesis and described in Chapter Three. The 
third objective can be achieved by combining the approximate eigenvalue 
assignment techniques to be described next with the exact eigenvalue 


assignment techniques. This can be done by using a crucial design tool 
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developed and employed in Chapter Three. Further information about 
this combined approach will be provided in the simulation study of 
Chapter Four. 

The output feedback algorithm developed in Chapter Three is 
based on some of the ideas involved in [31, 35, 45]. In addition to 
enabling the designer to assign more eigenvalues, it provides additional 
design freedom through its recursive nature which in turn may lead to 
more successful applications. Since it includes the algorithms of 
[31, 35, 45] as special cases, the discussion of these algorithms will 


be omitted from this chapter. 


Approximate Eigenvalue Assignment Techniques 


Approximate eigenvalue assignment methods do not guarantee 
the exact assignment of any of the system eigenvalues, and have been 
devised mainly to stabilize unstable systems. General information on 
the subject of output feedback stabilization can be found in [48, 49], 
and necessary and sufficient conditions for system stabilizability by 
constant output feedback are provided in [50, 51]. 

For approximate eigenvalue assignment, Jameson [45] has 
Suggested the use of a dyadic feedback controller, G = gk, to 
minimize the value of the closed-loop characteristic polynomial, r(AdA), 
when it is evaluated at the desired eigenvalue locations. This can be 
done by selecting g to fulfill the condition in (2.49) and determining 
a ki which will minimize the following objective function: 

r 2 


jeer ) Sir Oe te he (2.72) 
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Jameson [45] has generalized this objective function in order to weight 
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the different eigenvalues according to their importance and to minimize 
the weighted sum of the squares of the elements of ki, This results 
in the following objective function for the approximate assignment of 


all the system eigenvalues: 


5 Vilng: sacal ual ed“ (28735 
where L, represents the weighting on the pth eigenvalue and Be 
the weighting on the jo element of KI, Further information on 

a modification of this approach will be presented in Chapter Four. 

A similar approach has been suggested in [19, 20] and its application 
in a simulation study involving the control of a gas-turbine has been 
very briefly described in [52]. These methods have the disadvantage of 
not differentiating between a positive and negative deviation of an 
eigenvalue from its desired value which is important in the stabiliza- 
tion of continuous-time systems. This difficulty is due to the choice 
of the objective functional in the minimization problem. An advantage 
of the method is that it results in a rather simple expression for 


a The obvious disadvantage in its application is the lack of 


information about conditions which ensure satisfactory eigenvalue 
assignment and a satisfactory transient and steady state response of 
the system. 

Koenigsberg and Frederick's [49] and McBrinn and Roy's [53] 
approximate eigenvalue assignment techniques are gradient-search 
techniques based on the sensitivity of the system eigenvalues on the 
elements of the feedback gain matrix. Both papers [49, 53] attempt 


to shift the "poorly" located open-loop eigenvalues one-by~-one to 


better positions without shifting the others significantly. The 
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iterative procedure can be repeated as many times as desired although 
continual improvement of the results is not guaranteed. The sensitivity 
expressions used are due to Rosenbrock [54] and Reddy [55], who in turn 
have used the work of Fadeev and Fadeeva [56]. 

Still another eigenvalue sensitivity based method is that of 
McSparren and Etzweiler [57], whose method involves the change of some 
elements in the A matrix in order to approximately assign all of the 
system eigenvalues. They also suggest an iterative approach and provide 
some guidelines in the choice of state variables to be measured. A 
severe disadvantage of the method is that it is limited to only single- 
input systems which are completely controllable and observable. 

The successful application of these methods to the stabili- 

Za lon uoLeup to pee order systems have been reported [49, 53, 57], 
which is encouraging. 

Porter [58-60] has proposed the application of Liapunov's 
direct method in shifting the open-loop eigenvalues of a system to more 
stable locations thus stabilizing (or increasing the degree of stability) 
of a system. As noted by Johnson [61], Porter's method has extremely 
restricted applicability since it requires the existence of n indepen- 
dent measurements and controls, where n is the order of the system 


under consideration. 


2.10 Use of Observers and Compensators in Eigenvalue Assignment 


In systems where all of the state variables cannot be 
directly measured and the use of output feedback controllers does not 
provide adequate eigenvalue assignment, the designer has the option 


of using observers and compensators. This will ensure exact assignment 
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of all of the system eigenvalues provided that the system is completely 
controllable and observable. 

Reconstruction of all of the state variables by a reduced 
order Luenberger observer [62] and their subsequent use in a state 
feedback law ensures the assignment of the (2n-r) eigenvalues of the 
resulting closed-loop system. The separation property of the observers 
[62] allows the designer to separately specify the eigenvalues of the 
controlled system and the observer, and to design the observer and 
controller independently. 

The application of this approach to physical systems has been 
reported in the literature [63, 64]. Sturgeon and Loscutoff [63] have 
tried to control a double inverted pendulum both in simulation and 
experimental studies. Their model of this system is 6rR order with 
one control and three measurements, and their design approach consists 
of positioning the observer eigenvalues by the methods of [40] and 
those of the controlled system by the methods of [1, 5, 30]. Although 
simulation studies performed on an analog computer gave fairly good 
results, the experimental studies resulted in limit cycles [63]. 
Sturgeon and Loscutoff attribute this to modelling errors, the non- 
linearity of the real system and noise in the measurements. Their 
simulation studies clearly indicate that the observer rather than the 
controller is the source of the difficulties encountered. McMorran [64] 
in praising the Inverse Nyquist method [8] refers to Munro's Ph. D. 
thesis where an observer in conjunction with an eigenvalue assignment 
technique has been used. ‘He blames the controller for the extremely 
poor response of the controlled system during regulatory operation. 


Simulation and experimental observer studies at the University of 
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Alberta [65] gave similar results for an optimal controller. The poor 
responses were not due to the controller but to the observer which can 
be quite sensitive to unmeasured disturbances. 

The degree of the total system consisting of the controller 
and observer can be reduced if only some of the system eigenvalues have 
to be shifted. This can be achieved by the use of functional observers, 
which aim at reconstructing G x rather than x. This method has been 
studied in [66, 67], but not applied to an eigenvalue assignment problem. 
Although this reduction of the observer order is important because of 
the high order of chemical systems, it still does not solve the vital 
problems associated with unseen disturbances. Smith and Davison [68] 
have suggested the obervation of disturbances. This approach may be 
promising, but requires further development. 

The use of dual observers [62] is still another alternative 
for assigning the eigenvalues of the controlled system and the observer. 
In fact, the dual oberver design procedure devised by Murdoch [69] 
reduces the order of the observer to (Vorb, where we is the con- 


trollability index of the controlled sytem. 


Compensators 

Pearson and his co-workers [70-71] have suggested augmenting 
the dynamic system to be controlled by an ingenuously designed dynamic 
system such that the total system has the desired set of eigenvalues. 
The dynamical system to be added, i.e., the dynamic compensator, has 
the minimum degree of min (vocl, vy); where Wr and Vv, are 
the controllability and observability indices of the system to be 
controlled. The compensator design technique suggested in [70, 71] 


is a time domain method and provides for the assignment of all of the 
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eigenvalues of the augmented system. Thus, one does not have the option 
of specifying the eigenvalues of the controlled system and the com- 
pensator separately. Also, the designer cannot reduce the order of 

the compensator by requiring the assignment of only a subset of the 
system eigenvalues. 

Ahmari and Vacroux [72] have suggested a compensator design 
technique based on that of [70, 71] and the eigenvalue assignment 
technique of [31, 35]. They have shown that for every compensator 
order selected there exists a unique maximum number of assignable 
eigenvalues. Shaw [73], on the other hand, has noted that for the 
assignment of certain sets of eigenvalues of the total system, the 
compensator might have to employ unstable eigenvalues. This, in turn, 
drastically increases the sensitivity of the controlled system eigen- 
values to parameter variations in the system [73]. 

Chen [74] and Chen and Hsu [75] have proposed frequency 
domain compensators, whose minimal order is the same as in [70, 71]. 
The technique of [74, 75] allows the designer some freedom in con- 
figureating the compensator with the controlled system, Chen [74, 75] 
has considered two configurations, one of which allows independent 
specification of compensator and controlled system poles thus 
alleviating some of the disadvantages of the method of Pearson and 
his co-workers. But this extension may increase the order of the 
compensator. 

The methods of this section have apparently not been applied 
to physical systems. They possess some design Prcanecere which could 
be employed for design objectives in addition to eigenvalue assignment. 


The introduction of additional dynamics to already high order systems 
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is rather difficult to justify unless the added dynamics can be pur- 
posfully employed. 

A chief advantage of eigenvalue assignment techniques is that 
they are relatively easy to design and to implement. The use of 
observers and compensators in eigenvalue assignment tends to reduce 
this advantage considerably. 

An interesting interpretation of the methods described in 


this section can be found in [76]. 
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CHAPTER THREE 


A NEW ALGORITHM FOR EIGENVALUE ASSIGNMENT 


USING OUTPUT FEEDBACK 


3.1 Introduction 

In recent years the design of control systems which provide 
for arbitrary pole assignment using eatput feedback control has received 
considerable attention [1, 2, 6-9]. Several design methods have been 
reported which employ a control matrix which is a dyadic product of two 
vectors [1, 2, 8]. Typically, the number of eigenvalues which can be 
arbitrarily assigned is restricted to max (m,r) where m and r 
are the dimensions of the output and control vectors, respectively [2]. 
In this chapter, a new algorithm is presented which allows min [mtr-1,n] 


poles to be assigned subject to certain mild restrictions. The algorithm 
provides a simple, analytical expression for the output control matrix 


as the sum of two dyadic products. 


3.2 Formulation 
Consider the completely controllable and observable system 


described by: 


2) 


Ul 
Wa 
[x 


where: jx) ts Jan on, 1. sfategvector,, UG ts an, rox) (control vector, 


y is an mx 1 output vector, and A, B and C are constant matrices. 


Davison and Chatterjee [2] have given sufficient conditions for the 
arbitrarily close assignment of max (m,r) eigenvalues by the follow- 


ing linear, constant feedback law: 
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Top) 


Var C3ie2) 


A relation between the closed-loop and open-loop characteristic 
polynomials, r(A) and q(\), of the closed-loop system consisting of 


(3.1) and (3.2) has been derived by Hsu and Chen [3]: 


ECA) 
q(A) 


= det [I - CQI- a)" 8 


Wa 


dine (3.3) 


Assume that the feedback control matrix, 


lq 


iigteseor che 


dyadic form 
Giz goktros (3.4) 


The resolvent matrix, (Al - ye for any A #0 can be expressed 


as [4]: 
FO)s qtap + & POO * 4.0 + FP 
Cos ay Eten FADE te eA ec Sangha Aeceng eet 9 (3.5) 
af Ges: q(A) q(A) ; 
where Fi Fo» ae | can be obtained from Leverrier's algorithm 


(5) .- Combining .(3.3)—(3.5) “gives f1]* 
r() = 4Q) —k CFO) Bg . (3.6) 


The choice of a dyadic feedback matrix can be used to form 


two systems which are equivalent to the original system in (3.1) (3.2) 


and (3.4): 
System 1: Sea Xe De Ul 
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System 2: OSTA Ke baie 
(3.10) 
a 
Yee 
with ce =k ¢ (3.11) 
eae cay ae. C3. 12) 


It is obvious that Systems 1 and 2 have the same closed-loop 
eigenvalues as the original system. Also, Ding, Brasch and Pearson 
[6], and Davison and Wang [7] have proved the existence of g and 
ki which results in completely controllable and observable systems, 
respectively. 


For System 1 algebraic manipulation of (3.6) gives 


r(k) = (QA) - k* CQUe (3.13) 
and for System 2: 

r(A) = q(A) - 2. Bi SxUrer, (3.14) 
where Gre (Sihabeegh, 265 ane” Bt (3.15) 

g= eM” SUG ec, Lush” cel (3.16) 

1 a, a, a3 +s al 

Dine a, 4 a) 

Ue ne OE CO, yu Lita aba (Gi) 
0420 0 0 i 
n-1 n-2 n-3 fi 


|@ 
(pare | 
> 
w 


eS oe ee A (3.18) 


Saget 


qool-beseols omue sax ova S bas f amesahet said ie 


or 
noesess baa doesyvd cere ,celA  .medeye Lantgtro oe ea eta eg. 


bas 9g to sofetulye six bavorq svad [4] gnsW i nes t bor 
lary le hs : a 
,amsseye aidevrsedo: bas side lLloramoo Vee meee 


f: 
yore Plea +: * 
as SR, 
2svig (6.8) te aokteloginsm shendonis . oe ee 


ad 
~s 


hal an 
EL 56) aUe?2 Ma ~ Wp = Oy Sea 
we A) . 
One) . hee 
ve a 8 
Saale 
(21,&) {d “4 este 
oe Tent 
(OL.€) [> Ke BY ess at) 
; 
ten 
ip. 
Cog 
vas . Bb 
eh 
oT : 
L f 
T "8 
(61.2) pL ttre a 


54 


and aj» a a are the coefficients of qQA) , which can be 


22 © 0 2 9 n-1 
determined from Leverrier's algorithm [5]. Thus (3.13) and (3.14) give 


the relations between the open-loop and closed-loop eigenvalues and 


the controllability and observability matrices, Q and 


ihe) 
rv 
Oo 
Fh 


Systems 1 and 2. 

Davison and Wang [7] have shown that almost any output feed- 
back control matrix will generate a closed-loop system with a set of 
eigenvalues disjoint from that of any corresponding completely con- 
trollable and observable open-loop system. Thus, in what follows, it 
will be assumed that any open-loop system (or equivalent system) either 
possesses distinct nonzero eigenvalues or has been converted into one 


with distinct nonzero eigenvalues. 


3.3 Algorithm 


In the derivation of the eigenvalue assignment algorithm, the 
two casés of m>r, J‘and” % >\miwitl be’ considered separately. \iIn 


both cases the control law is expressed as 


usu, tu, (3.19) 
h = as (3.20) 
where Dees cag : 
Pet comeciee (3.21) 
Mo * £9222 ° 


Gaseck: ome © 

The algorithm involves three steps: assignment of m 
distinct eigenvalues, "protection" of (r-1) of the eigenvalues (i.e., 
making them invariant to subsequent feedback control), and assignment 
of m additional distinct eigenvalues. Thus the total number of 


distinct eigenvalues assigned will be mtr-l1 (or n if mtr-l >n). 
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Step 1 (Assignment of m_ eigenvalues) 


Choose a gy vector such that the resulting equivalent system 


of (3.7) becomes: 


(3222) 


ie ‘, Sone 
where A, =r; by = Bg) and u_ = ky y. Evaluate (3.13) at m 


r r 


desired eigenvalues, i dpa om eae 


dl? and impose the conditions 


that 


1, O4)) = 6, O42) =... = 1 Og) = 0 (3.23) 


where r is the closed-loop characteristic polynomial for the system 


in (3.22). Then k, can be determined from (3.24) below which is 


1 


similar to Davison's algorithm [8]: 
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that dr r belong to the set of eigenvalues of the closed- 


div’ da? i toda 
loop system matrix (A, + B g) ki C)e es ie=invertibiiaty ‘of (C Q. UE) 


is determined by the choice of the desired eigenvalues. 


Step 2. (Protection of (r-l) eigenvalues) 


Consider the closed-loop system resulting from Step 1 as a 


new open-loop system given by: 


ae ao 2eeE BS 
(B27) 
yO eee ee 
where 4, = A. + 3B Kk? C 
22.) #1,)= 41-1 =: 


Output feedback affects all of the sytem eigenvalues which 
correspond to both controllable and observable modes [7], while none 
of the eigenvalues belonging to either an uncontrollable or unobservable 
mode can be affected by any output feedback controller [1], [9]. Thus, 
a subset of (r-l) eigenvalues of (3.27) can be "protected" by introduc- 


ing an equivalent single-input system: 


x + db. 


2 Zee 


(3.28) 


with b, = B and u, = ic? y, in which the modes corresponding to 
these (r-l) eigenvalues have been made uncontrollable by an appropri- 
ate choice of Bo: However, in order to assign m additional eigen- 
values in Step 3, it is also necessary for the system in (3.28) to con- 
tain at least m controllable and observable modes. (Observability 
of the sytem in (3.27) is not affected by going to (3.28).) These con- 


ditions can be achieved by choosing & such that 
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and V (3.30) 
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where hy is an arbitrary vector with at least m nonzero elements 
i ; 

and V5 contains as its rows, (r-l) of the left eigenvectors of A, 

corresponding to (r-1) of the m eigenvalues which were assigned in 

2 

2 


Step 1, while V, has as its (n-r+l) rows, the remaining left eigen- 


vectors of A: 
Equation (3.29) can always be satisfied, since it requires 

finding an r-vector, Bos which is orthogonal to (r-1) r-vectors. 

bop tack, If rank (V; B) < r-l, then the vector, cD) which satisfies 

(3.29) is not unique. The condition in (3.30) requires that the choice 

of B must be nonorthogonal to at least m vectors, which are a sub- 


set of the rows of v2 B. (The number of rows of v B which fulfill 


2 
this condition is equal to the degree of the relative minimal polynomial 


of b, for the matrix A MO1.) 


2 2 
Remark: In situations where a particular choice of cD) 

satisfies (3.29) but not (3.30), several alternatives are available. 

One could return to Step 1 and choose a different &)> which will, 

in general, change the (n-m) umassigned eigenvalues of (A, wee) 8 kG) 

and consequently result in a new (Ws B) matrix, or one could perturb 

the desired eigenvalues slightly and repeat Step 1 hoping that the 

resulting system will allow (3.29) and (3.30) to be satisfied. Alter- 


natively, one could add another feedback controller with feedback gain 
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to (3.27) thus affecting its (n-rtl) eigenvalues and all of the system 


eigenvectors. 


Step 3. (Assignment of another m_ eigenvalues) 


Consider the equivalent system resulting from Step 2 as a new 
open-loop system given by (3.28). Application of Step 1 to this system 


will yield 
T -1 
ky = qd) (€ Q, U, E,) (3232) 


where Q,> U,> E. and 4d, pertain to the system in (3.28) and are 
analogous to the quantities defined in (3.15), (3.17), (3.25) and (3.26), 
respectively. It follows that rank (U.) =n, rank (E,) =m and 

nN > rank (Q.) >m. The last inequality results from the sufficient 
condition set for successful application of Step 2, and the fact that 

the rank of Q, is equal to the number of controllable modes of the 


system in (3.28) [10]. Thus the condition 


rank (C Q, U =m (333) 


is sufficient for the successful application of both Steps 2 and 3. 


Case 2: m<r 

The algorithm in this case is the dual of that described in 
Case 1. It involves: assignment of r eigenvalues, protection of 
(m-1) of these, and assignment of r additional eigenvalues. Thus, 
the total number of eigenvalues assigned is (mtr-1) provided that 
(ortr=1) < mn. 

Here the eigenvalue assignment is based on the use of the 
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Step 2 which makes (m-1) of the r modes (whose eigenvalues have been 
assigned in Step 1) unobservable while leaving at least r of the 
(n-mtl) remaining modes observable. In analogy with Case 1 the final 


control law becomes: 


u=g, ki + g& k, (3.34) 
where 
ley =i eT. 
ete (Ey U, 3) B) a7 (3.35) 
OT TT —-f T 
8 ~ (E U, 8 B) qd» (3.36) 
and 1° S»> UL and U, are defined in analogy with (3.16) and (3.17) 
at each stage. k must fulfill the following conditions: 
Ay i 
ko W = 2 33) 
ik Z 
ky OW) = £ (3.38) 
where the m-vector £, contains at least r nonzero elements. WS 
has as its columns (m-1) of the right eigenvectors of A, and ws 


has as its columns the remaining (n-m+l) right eigenvectors. 
Cases 1 and 2 then form a constructive proof for the follow- 


ing proposition: 


Proposition 


If the system in (3.1) is completely controllable and 
observable with rank C =m and rank B = r, then there exists a 
constant output feedback matrix G such that min (mtr-l,n) eigen- 
values of the closed-loop system can be assigned to arbitrary, distinct 


locations (with complex values chosen as conjugate pairs) provided that 
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3.4 Example 
Consider the following third order system with two control 


variables and two output variables: 


Xp -3 2 0) xy 0 a us 
X, = 4-5 3 ap et |i (0 u, 
X, OF 0) 3 X, 0 i 
(3.41) 
it f 1 0 0 x 
Yo 0 if 0 Xo 
ee 
Its eigenvalues are given by 
A, = -1 1, = -3 DE pet hi (3.42) 


and it is desired to find a control law, u = 8) ky SP & Lee such 


that the closed-loop eigenvalues become: 


dad = -10 d 42 = ~9 43 = -8 ., (32.435) 
Choose 
ut 
& ms (3.44) 
iL 
then k. can be determined from (3.24) as 
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ky = [-8.27 O.09}) (3145) 
cE 
The resulting system equation, A, = A, eB gy kK C, has as its 
eigenvalues: 
X,, = -10.00 Ay = -5.00 OP cet al ho Maas 


dl d2 d3 


Now, choose a B which satisfies (3.29) and (3.30) as 


DETTE 


& ° (3.46) 
1.000 


k, can be calculated from (3.32) as 


ji 
De 
k, = [-96.73 34.58] (3 41) 
which yields the final closed-loop system matrix, A, = A, + B 8 k, C» 


whose eigenvalues are 


ad = -10.00 42 = ~9.00 43 = -8.00 . (3.48) 


It is worth noting that if previously reported design pro- 
cedures [1, 2, 8] were employed, only max (m,r) = 2 poles could be 


arbitarily assigned. 


3.5 Conclusions 

A new algorithm has been presented to achieve pole assignment 
using a proportional output feedback control law. Provided that certain 
mild conditions are satisfied, the numbér of poles that can be assigned 
to arbitrary, distinct locations is min [mtr-l,n] where m, r and n 


are the dimensions of the output, control and state vectors, respectively. 
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This number is a significant improvement over the max (m,r) poles 
which can be assigned using previously reported algorithms [1, 2, 8]. 
The algorithm consists of a two stage design procedure where 
some of the closed-loop eigenvalues which are shifted during the first 
stage are then "protected" by transforming the system into an equivalent 
single input (or single output) system. The equivalent system is chosen 
so that it satisfies the conditions of "ideal control" (or "ideal 
measurement") of Takahashi et al [12] since some of the modes have 
been made uncontrollable (or unobservable). Application of the second 
design stage allows assignment of additional poles without affecting 
the "protected" poles. 
A numerical example demonstrates that the algorithm is quite 
simple and results in an analytical expression for the output control 


matrix. 
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CHAPTER FOUR 


SIMULATION RESULTS FOR A DOUBLE 


EFFECT EVAPORATOR 


4.1 Introduction 


Some of the modal control methods described in Chapter Two 
and the eigenvalue assignment algorithm presented in Chapter Three 
have been applied to the fifth and third order state space models of 
the pilot scale, double effect evaporator in the Department of Chemical 
Engineering at the University of Alberta. 

In this chapter the author has attempted to describe the 
design options available in the various design methods and to denm- 
onstrate their usefulness in fulfilling two design objectives in 
addition to specification of the closed-loop eigenvalues. These 
design objectives have been minimization of offset in the most 
important process variable, the product concentration, and some 
control over the magnitude of the gain elements in the controller 
matrices. These objectives have been deemed important, since inte~ 
gral action was not incorporated into the controllers used and 
excessively high gains would lead to difficulties in the future 
experimental evaluation of the controllers designed in this chapter. 

The simplicity of the resulting controllers and the modest 
amounts of computer time required to design them have enabled the 
performance of a great number of controllers to be evaluated; only 


a few of these will be described here. 


63 


Vy iF e que 


“agua A got pe ter. 


at 


serit raiqen) al beineesxqd mAasivogis parse 


to elehom aos02 sisie byo bibkds bas dahee ads 03 bet ;  ased 
: 25 


isolmenm 20 tItaqgoll sft at tosatogsys as336 sldvob .olnee 2 


I -. 7 a 
.BIvedIA Yo Wievevinw ots 48 ge sine 


—_ 4 _ 
= 


edi adizsash oF bsjqmesta aac yodsua siz sed qardo elds ’ 

~meb o3 bas sbortsem agitesb avolusy sd3 al aldsiisvs — enotsac 
- ; a ic a = 

al ssvisgsetdo agteab ows gallitiiu? at essniuisey sisi os bg 

ae 


seotT ,.asulsvesgke qool-bsaols sit to nolssotttooge: oF ac ove 


Jeom Sci Gk Seetto to notvesiImiata ceed ovad evtiostdo ngtes. 


icy Tis 


omj08 hos 5 


,aotseyinsonos towborg sf «olds tev e280 


4 


t9iforsnos of3 mt asnemots tiiag o3 to obuttngem ers 


~etnt sonts ,jnstdxoqgml bemesh ased evar mnie fi 
te 

bos bseu etelior2nos si2 oat heseseqreaah; — 

stujut edi ot estsiuolttitb o3 bsel vied 

| - 


@ 


-tesqendo ald3 ot benghesh erslfoxinos an to sotjaul BN eee 9) 
aa Liat 
‘ -~ 49 rae 
tasbom sfi3 bas eteliortn0. gatiivess ed to bare “oat 


es 
sin’ ents a duqmos to. 


4 


oid bsidans eved meds ngteeb ot bas 


- a 
de 3ss1g 8 to soneartot 


ae ar ed 
ediupeeb ¢ £ iw aod: 


Vino thesauieve ed o3 ¢ ‘eonisieabild ea 


ee Soe 

aT. TF 

a oe = 
a 


64 


4.2 Physical Description and Mathematical Model of the Evaporator 

Some of the control algorithms described and developed in 
Chapters Two and Three have been applied to the pilot plant scale, 
double effect evaporator schematically represented in Figure 4.1. 

The first effect is a natural circulation calendria type 
unit which is fed with 5 1b./min. of 3.2 percent by weight of 
triethylene glycol. The feed is heated by 2 1b./min. of fresh steam. 
The second effect is an externally forced - circulation long tube 
vertical unit, which is heated by the first effect vapor. It:pro- 
duces about 1.5 1b./min. of ten percent triethylene glycol. The 
second effect vapor is totally condensed. Tight pressure control 
maintains the necessary pressure differential between the effects. 

The systematic modeling of the evaporator under consider- 
ation and of similar units has been extensively studied at the 
University of Alberta [1-4]. Among the models developed for the double 
effect evaporator, a tenth order continuous-time nonlinear model 
derived from material and energy balances has been found to represent 
the experimental open-loop system behavior best [3]. 

The application of state space control methods devised for 
linear systems has necessitated the linearization of the nonlinear 
model. Open-loop experimental studies have indicated that this 
linearization introduces a significant modeling error unless the 
process variables stay roughly within a ten percent range of the 
steady state values around which linearization has been performed [3]. 
Regulatory control aims at holding the process variables close to 
their steady state values; thus, modeling errors introduced through 


linearization are expected to be smaller for the closed-loop system 
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Both the design and implementation of multivariable control | 


techniques can be carried out more easily and more economically with 


lower order models provided that these models predict the system 


behavior satisfactorily. The tenth order linearized evaporator model 
has been reduced to fifth and third order models through physical and 
modal considerations by Newell and Wilson [2-4]. Open-loop simulation 


and experimental studies have suggested that the use of reduced order 


models for multivariable control studies rather than the tenth order 
linear model is jas tl eiab Ve. This is due to the fact that the neg- 
lected modes of the higher order model are considerably faster than 
the retained ones. 

The fifth and third order disorecent ime evaporator models 
used in this chapter have been derived by Wilson [4]. They are of 


the form: 


Sine) Lig =g extn) teas St O° (nt) C45))) 


y (nde 


Wa 


where T denotes the time base or discretization interval. Wilson 


“(nT)” « (4.2) 


has applied Marshall's model reduction technique [5] to the tenth order 


linearized continuous-time model to obtain the fifth and third order 
continuous-time models, which then have been discretized wich a time 
base of T= 64 seconds. This discretization is desirable for 
digital simulation and direct digital control applications. 


Theo, 45 0° and © (matrices tor the fifth and third 


order discrete evaporator models are presented in Tables 4.1 and 4.2, 
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TABLE 4.3 


Description of the Evaporator Variables 


plate Vector, x 

x = Iw1, Cl, H1, W2, @2] 

Wl First effect holdup 

Cl First effect concentration 
Hl First effect enthalpy 

W2 Second effect holdup 


C2 Second effect concentration 


Control, Vector, ».v 

uw = [s, Bl, B2] 

S Steam flow 

Bl First effect bottoms flow 


B2 Second effect bottoms flow 


Disturbance Vector, d 
As 

dontilEs,GE, HE] 

F° Feed flow 

CF Feed concentration 


HF Feed enthalpy 


Output Vector, y 
y. =, [w1, w2,, c2] 


Normal Steady State Value 


4G. Splibis 

4.59% glycol 
189.2 BIU/1b. 
41.5 -1b. 


10.11% glycol 


2.0. 1b..., min « 
3.485 1b./min. 


5845 Lb ./min. 


er bas) Mat. 
3.2% glycol 


156.9 BIU/1b. 
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respectively. The state, output, control, and disturbance vectors are 
defined in Table 4.3 and are expressed in normalized perturbation form, 


a 


_ W1-Wiss 


q 
we Wlss 


(ss = steady state) , (4.3) 


Their normal steady state values are also presented in Table 4.3. 
These same physical variables are involved in the third order model 
except that only W1', W2' and C2' have been retained as the state 
variables. 

The accumulated experience on the experimental behavior 
of the evaporator over several years suggests that its closed-loop 
behavior deteriorates when the gain elements in the feedback matrices 
are excessively large. This has been attributed to neglecting noise 
effects, nonlinearities and time delays. The recommended upper limit 


on the value of gain elements is 100.0. 


4.3 Modal Analysis of the Discrete Fifth Order Evaporator Model 

Some of the design methods discussed in Section 2.6 will be 
applied here to the fifth and third order evaporator models. The 
eigenproperties of these models are in Table 4.4. The arrangement of 


the matrix element is such that 
V 0) W = A ‘ (4.4) 


Table 4.5 present the mode controllability and mode observability 
matrices for the fifth order model. 
Consideration of Tables 4.4, 4.5 reveals that: 


1) The open-loop system has two unstable modes (corresponding 
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to eigenvalues of 1.0). 

2) The unstable modes possess repeated eigenvalues with 
distinct eigenvectors. 

3) All of the system modes are both controllable and 
observable. 

4) The system matrices are almost diagonal. 

5) The right and left eigenvectors of the system are fairly 
well distributed in the state space. 

The instability of the open-loop evaporator is due to the 
integrating nature of the two liquid levels involved. This type of 
instability is rather common in chemical engineering systems, and 
thus of considerable interest. Any satisfactory control systems 
design must be able to stabilize these unstable modes. 

The presence of a pair of repeated eigenvalues with distinct 
eigenvectors implies that a unity rank feedback matrix cannot stabilize 
the system. This follows from well-known controllability consider- 
ations. In order to shift both of the unstable eigenvalues, the 
designer must either increase the rank of the feedback matrix or add 
an "insignificant" feedback control loop [6] to generate a new open- 
loop system with distinct eigenvalues. Alternatively, one could per- 
turb the elements of $ slightly so that @¢ has distinct eigenvalues. 

The approximately diagonal nature of the system matrix may 
imply that the eigenvectors form almost orthogonal sets, i.e., the 
right eigenvectors are almost orthogonal to each other and so are 
the left eigenvectors. Consideration of Table 4.4 indicates that 
this is not quite true. The reason for this is the high numerical 


sensitivity of the eigenvectors to the elements of the system matrix. 
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Thus, by considering the system matrix alone one cannot decide on the 
relative influence of each mode on each physical variable. But, with 
the eigenvectors well distributed in the state space it is not dif- 
ficult to see that: 

1) Wl is most influenced by the first mode, i.e., the mode 
whose right eigenvector of the first colum of W in Table 4.4. 

2) C2 is most influenced by the second mode. 

3) W2 is most influenced by the third mode. 

It is also apparent that one cannot control the third mode without 
affecting the first mode considerably. Fortunately, this and similar 
interactions among the system modes have not been so severe as to pro- 
duce unstable modes in any of the simulation studies to be described 
in the subsequent sections. 

Next consider the open-loop response of the system to separ- 
ate twenty percent step changes in feedflow, F, feed concentration, 
CF, and feed temperature, TF shown in Figure 4.2. The unstable 
response of the two holdups, Wl and W2, is apparent. The dis- 
turbance caused in the product concentration, C2, the main variable 
of interest is also significant. Thus, for any control strategy to 
succeed the first, the second, and the third modes must be controlled. 
This is also the reason for keeping these modes in the third order 
model. 

The small dimensions of the system under consideration makes 
the configuration of the controller quite simple. The three measure- 
ments to be taken are obvious, and so are the three physical variables 
to manipulate in order to control W1, W2 and C2 effectively. 


The fact that the eigenvectors are well distributed in the 
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State space has two implications on the performance of the modal con- 
trollers to be designed: 

1) The modes will not be highly interactive. 

2) The magnitude of the controller gains required for a 
certain eigenvalue shift using modal control are expected not to be 
significantly greater than that required for any eigenvalue assign- 
ment technique. 

The basis for the last argument can be found in McFarlane's 


interpretation of modal control [7]. 


4.4 Application of Rosenbrock's Approximate Modal Control Method [8] 
to the Fifth Order Evaporator Model 

Rosenbrock's modal control technique has been simulated in 
this section. This technique has the desirable feature of providing some 
a priori indication of the success to be expected from its application 
to a specific physical system. This a priori evaluation for the 
evaporator model has been presented in Section 4.3. It is worth noting 
here again that the open-loop system eigenvectors are directionally 
well distributed in the state space and the eigenvalues corresponding 
to the dominant modes of the system are fairly well separated from the 
rest of the eigenvalues. In fact, the good distribution of the eigen- 
vectors implies that each state variable is affected most by one 
particular mode which is different for each state variable. 

The information about how well the eigenvalues of the control- 
led modes will be assigned can be obtained by applying Takahashi et al's 
analysis [9] described in Section 2.7. Application of their analysis 


[9] using the mode controllability matrix, H, and the mode observability 
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matrix, F, given in Table 4.5 suggests that the approximate modal 
controllers to be designed can be expected to assign two of the system 
eigenvalues quite well. Also the eigenvalues belonging to the wncon- 
trolled modes will not be moved significantly from their open-loop 
locations due to modal interactions. Specifically, sede vlan 
by the zero elements of the last two rows of the H matrix and the 
Structure of the modal-domain closed-loop system matrix of (2.65). It 
should be noted though that the above arguments are only approximate, 
and unless either (2.62) or (2.63) holds, the exact assignment of any 
eigenvalue cannot be guaranteed. 

The availability of three measurements and three controls 
implies that three modes can be approximately controlled [8, 9]. But, 
as indicated by the application of Takahashi et al's analysis [9] to 
the evaporator model, it will be possible to assign at most two of the 
three eigenvalues of the poneEo Ted modes 

After having decided on which modes to control as described 
in Section 4.4. The following design options are available and will 
be investigated in the simulation study: 

1) The order in which the eigenvectors corresponding to the 
controlled modes appear in Wo and ve which are used in the cal- 
culation of the ‘controller,matrix jG via CZaG ee 

2) The specification of the desired eigenvalues, Aa? 

3) Pairing of the open-loop eigenvalues with the desired 
closed-loop eigenvalues. 

In the remainder of this section two major objectives of the 
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To demonstrate the usefulness of the above design options as 
design tools and to correlate the system closed-loop response, the eigen- 
value locations, the eigenvector directions, and the characteristics of 
the controller matrices. 

It should be noted that in cases where ideal modal controllers 
can be designed, only the second design option (above) has an effect on 
the resulting modal controller and thus on the performance of the 
resulting closed-loop system. But, in all other cases these design 
options, in general, can be used to obtain different controller matrices, 


thus different closed-loop eigenvalues, eigenvectors, and transient 


responses, 


Simulation Study 

In the simulation studies to be described in this and the 
following sections, disturbances in feed flow and feed concentration 
presented comparable challenges to the controllers, while disturbances 
in feed temperature were less severe. Thus, most of the simulation 
studies of this chapter will involve feed flow and feed concentration 
disturbances only. Also, the performance of all of the controllers 
will be evaluated by considering 20% step changes in feed flow and 
feed concentration beginning at time t = 0. This provides a common 
basis for the various comparisons that are made. 

The eigenvalue and eigenvector calculations in this thesis 
have been performed on an IBM 360/67 computer. For this purpose the 
double precision version of the University of Alberta Computing 
Center subroutine CS201, based on Wilkinson's QR algorithm [10], 


has been used. The determination of the controller matrices and 
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transient responses have been carried out on an IBM 1800 computer. 

Table 4.8 reveals that the three last closed-loop eigenvalues 
resulting from all of the runs are fairly similar. Furthermore, inspec- 
tion of Table 4.10 indicates that the product concentration C2, which is 
the most important controlled variable, is affected most by the two 
modes corresponding to the complex conjugate eigenvalues. This is one 
explanation for the small differences observed between the C2 res- 
ponses resulting from the various runs. 

The location of the closed-loop complex eigenvalues is fairly 
close to the origin of the unit circle in the z-plane. This accounts 
for the well-damped behavior of the transient responses observed in 


every run. 


bitect of Elgenvector, Order 


The effect of the order in which the controlled mode eigen- 
vectors appear in the eigenvector matrices is considered in Runs 1-6. 
The details of these studies are illustrated in Tables 4.6 - 4.10. 
Figures 4.3 and 4.4 compare the best and worse responses that were 
obtained. The difference in the closed-loop eigenvalue locations in 
Table 4.8 is admittedly not great, although the controller gains for 
Runs 1-6 differ significantly (cf. Table 4.9). Two of the eigenvalues 
are fairly well assigned, as would be expected from the earlier dis- 
cussion. 

The effect of the order in which the uncontrolled mode eigen- 
vectors appear in the eigenvector matrices is investigated in Runs 7 
and 8. The resulting closed-loop eigenvalues and controller matrices 
were identical. This is entirely expectable since the uncontrolled 


mode eigenvectors do not have any effect on the controller matrix 
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designed. 


Desired Eigenvalues 

In choosing real desired closed-loop eigenvalues two points 
must be considered: the absolute and the relative magnitude of the 
eigenvalues. Runs 9-12 consider the first point. It is apparent from 
Figures 4.5 - 4.7 that the transient response improves with decreasing 
magnitude of the eigenvalues. Run 12 is not shown since its transient. 
responses were essentially identical to aaa ead Run 11. Consideration 
of Table 4.7 and Figures 4.5 - 4.7 indicates that this change is not 
directly proportional to the decrease in the magnitude of the desired 
eigenvalues. Also, comparison of Runs 1] and 12 clearly indicates 
that the amount of improvement to be expected by decreasing the value 
of the eigenvalues reaches a limit. Comparison of the corresponding 
closed-loop eigenvalues and eigenvectors indicates that there is only 
a small change in the eigenvector directions and also in the complex 
conjugate eine teiece eigenvalues which belong to the closed-loop modes 
most Soren C2. This suggests the importance of the closed-loop 
modes rather than cogs ate closed-loop eigenvalues, 

Similar observations can be made about the values of the gain 
elements in the corresponding controller matrices: they appear to tend 
toward limiting values as the desired eigenvalues become smaller. 
Another observation Orieone practical value, especially in ths par- 
ticular case study, is that the ner case in the eaves of the gain 
elements has not caused ary nonminimum phase behavior, which might have 
occurred according to MacFarlane [7]. Similarly, the largest gain 


elements are considerably smaller than the upper limit of 100.0 set 


in Section 4.2. 
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The second point, concerning the relative magnitude of the 
desired eigenvalues, has been illustrated in Runs 10, 13-15. Although 
the sum of the three desired eigenvalues is comparable the resulting 
controller matrices and the closed-loop responses in Figures 4.8 and 
4.9 are significantly different. Comparison of Runs 7, 12 and 17 . 
indicates that the relative magnitude of the eigenvalues has very 
small influence on the values of the gain elements in cases where all 
of the desired eigenvalues have small magnitudes. The runs considered 
in this paragraph demonstrates how important this design option can be 
in influencing the magnitude of the various gain elements, particularly 


in cases where some of the gain elements become undesirably large. 


Pairing of Open-loop and Desired Closed-loop Eigenvalues 


The last design option, pairing of the desired closed-loop 
eigenvalues with the open-loop eigenvalues, has: been considered in 
Runs 15, 18 and 19. Run 19, where the smallest desired eigenvalue was 
paired with the eigenvalue of the mode influencing both Wl and W2, 
resulted in the best controller while the controllers in Runs 15 and 18 
were quite similar as demonstrated in Figures 4.10 and 4.11. Thus this 
design option should not be ignored although it is not obvious a priori 
which pairing will result in the best combination. 

A closer look at the controller matrices given in Table 4.10 
suggests that all of them are quite similar. It is very difficult to 
interpret multivariable control laws by considering the individual 
gain elements in the controller matrices. But, extensive experience 
with the evaporator system and the multi-loop controllers designed in 
previous investigations [11, 12] indicates that each gain element in the 


controller matrix has’ a physically reasonable sign and magnitude with 
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the exception of the small elements in the (1,2) and (2,2) positions 
which tend to change sign from run to run. 

Since driving the modal activations to zero should eventually 
result in driving the state variables to their steady state values, the 
behavior of the manipulated variables will be physically meaningful as 
long as the modes are not highly interactive. This has beeen observed 


in all of the simulation studies of this section. 


Summary 


The following observations were made concerning the application 
of Rosenbrock's approximate modal control method [8] to the fifth order 
evaporator model: 

1) The product concentration, C2, which is the variable of 
greatest interest, showed a small offset and was well-behaved. The 
other state variables were also satisfactorily controlled. 

2) The manipulated variables, S, Bl and B2, showed in 
every run the physically expected behavior and never exceeded the 
physical limits of the actual evaporator. 

3) Two of the system eigenvalues could be fairly well 
assigned. Two of the other eigenvalues degenerated into a complex 
conjugate pair which changed very little from run to run. 

4) As the desired eigenvalues approached the limit zero 
the closed-loop eigenvectors approached limiting directions. Similarly, 
the gain elements in the controller matrices approached limiting values, 
which were reasonably small. 

9) The closed-loop dominant modes were those pertaining to 
the complex conjugate eigenvalues. Since these eigenvalues could not be 


directly affected the corresponding responses were similar in every run as 
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noted in point l. 

10) The choice of desired closed-loop eigenvalue locations, 
ordering of the controlled mode eigenvectors in the eigenvector matrices, 
and the pairing of open-loop and closed-loop eigenvalues proved to be 
useful design options. 

It is very difficult, if not impossible, to generalize the 
details of the evaporator results to other systems. But it is possible 
to draw some general conclusions: 

1) In systems with well-distributed open-loop eigenvectors 
and well separated eigenvalues, Rosenbrock's approximate modal con- 
troller [8] is expected to give satisfactory results. 

2) If one can also manage to design a control law with 
distinctly dominant closed-loop modes with desirable eigenvalues, the 
performance of the closed-loop system is expected to be very satisfac- 
tory unless the remaining closed-loop eigenvalues possess very undesir- 


able locations. 


4.5 Application of Takahashi et al's Approximate Modal Control Method 
[9] to the Fifth Order Evaporator Model 


Takahashi et al's approximate modal control method [9], which 
has been described in Section 2.7, lacks a theoretical basis for pre- 
dicting its suitability for constrolling a specific system. Thus it 
seems rather pointless to attempt a detailed study for a specific sytem. 
Consequently, only representative results were sought in this investi- 
gation. 

All of the design freedoms described and evaluated in Section 


4.4 exist for Takahashi's method, too. But, here only the third 
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design option has been evaluated. Experience with the various options 
in Section 4.4 has suggested that this option has the greatest 
influence on the response characteristics of the closed-loop system. 
The order in which the eigenvectors appear in the eigenvector matrices 
has been selected as the order used in Run l. 

Takahashi et al's method [9] allows the option of assigning 
all of the system eigenvalues. This was attempted in Runs 20 and 21. 
Table 4.13 indicates that although the eigenvalues have generally been 
shifted in the appropriate direction the eigenvalue assignement is 
very poor. The response characteristics in Figures 4.12 and 4.13 are 
also worse than any of the ones obtained in Section 4.4. In Runs 22 
and 23 two of the desired closed-loop eigenvalues were specified to be 
the open-loop eigenvalues of 0.1921 and 0.438 (i.e., only three of the 
open-loop eigenvalues were to be changed). 

Table 4.12 and Figures 4.12 and 4.13 clearly indicate that 
although five eigenvalues are shifted this strategy has significantly 
improved the response characteristics since the response times and 
offsets are reduced. 

The results in Table 4.13 indicate that it is very difficult 
to see a trend in how well the eigenvalues have been assigned. Sim- 
ilarly, it has not been possible to observe a trend in the closed-loop 
eigenvectors, which are therefore not reported. Consideration of 
Table 4.14, on the other hand, is very informative: Runs 20-23, 
although performed for different specifications of the desired eigen- 
values, have resulted in very similar controller matrices. This may 
have a mathematical explanation since these controller matrices rep- 


resent the least squares solution of a set of linear algebraic equations. 
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Apparently these equations are not very sensitive to the desired eigen- 
values used in Table 4.11. It is possible that use of the several 
design freedoms described in Section 4.4 can significantly affect the 
value of the gain elements and could help to improve the poorer response 
characteristics observed. 

A closer look at the controller matrices pertaining to the 
first four runs of this section indicates that gain element (1,3) 
an SG relating product concentration," C2, and’steam flow rate, Ss, 
is much smaller than those of the previous section. In Run 24 this 
gain element has been assigned a value of -10.00. The drastic improve- 
ment in the response characteristics (cf., Figures 4.12 and 4.13) 


indicates that intuition may be helpful, even in multivariable control. 


4.6 Application of Model Reduction and Ideal Modal Control to the 
Fifth Order Evaporator Model 

In Section 2.9 it was noted that model reduction offers an 
alternative to the application of approximate modal control techniques. 
Specifically, the designer can generate a reduced order system model 
with the same number of state variables as control and output variables 
which then can be used to design an ideal modal controller. The 
application of the ideal modal controller to the original model will 
form the real test for the justification of the approach used. Many 
model reduction techniques have been proposed in the literature. But, 
since the resulting model will be employed to design an ideal modal 
controller, Marshall's [5] model reduction technique based on the modal 
characteristics of the system has been deemed the most suitable. 


Reduction of the fifth order evaporator model to the third 
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order model described in Section 4.2 generates a system with 3 states, 
3 outputs, and 3 controls. Thus, the application of the ideal modal 
control technique described in Section 2.3 is possible. 

Two simulation studies were performed in which the desired 


eigenvalue sets were: 


{0.0000, 0.0000, 0.0000} (4.5) 
and 


10.0006" (0.0005, "O20004)) . (4.6) 


The resulting third order closed-loop matrices were exactly assigned 
these eigenvalues and the closed-loop eigenvectors were identical to 


the open-loop ones. Similarly, the resulting controller matrices 


were: 
edie) .0002 -13 .63 
Gt geo. 9.0000 5.465 (4.7) 
i223 24.61 24.14 
and 
1,039  =70491 -13.63 
Gea) eat co BeOS ot 509 ; (4.8) 
LZ. 7 24.79 24.13 


At this point one might note the similarity of these controller matrices 
to each other and especially to those obtained in Section 4.4 for small 
eigenvalues. 

Two tests were applied to evaluate the success of this 
approach: 


1) Determination of the eigenvalues of the fifth order closed- 
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loop system matrix resulting from the use of the controller matrices in 
(4.7) and (4.8). 

2) Determination of the closed-loop response of the fifth 
order model obtained by the controller matrices of (4.7) and (4.8). 

The results of these tests were very satisfactory. The 


eigenvalues of the fifth order closed-loop system matrices were: 


{O2001% O00 1s On 5 7020. 50717, 0.898} (4.9) 
and | 


{0.002, -0:005, (0.57720.561i,. 0.898) . (4.10) 


They are by no means identical to the eigenvalues of the third order 
closed-loop system matrices, but they are certainly much better than 

the open-loop system eigenvalues. The closed-loop system response 

to F, CF, and TF disturbances was almost identical using these two 
controllers. As depicted by Figures 4.14, 4.15 and 4.16 the closed-loop 
responses obtained in this section are better than any of those obtained 
in Sections 4.4 and 4.5. The offsets resulting from the two studies of 
this section were -0.002, 0.007, and 0.004 for F, CF, and MIF 
disturbances, respectively. 

The success of this combination consisting of a modal approach 
to model reduction and ideal modal control can by no means be generalized 
to other applications. Specifically, one cannot claim that this 
approach is superior to the approximate modal control techniques of 
[8, 9]. But, it is possible to say that, whenever the open-loop system 
possesses well separated eigenvalues and well distributed eigenvectors, 
the chances of this approach providing satisfactory results are quite 


good. 
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4.7 Application of the New Eigenvalue Assignment Method to the 
Fifth Order Evaporator Model 

In this section the eigenvalue assignment algorithm developed 
in Chapter Three has been applied to the fifth order evaporator model. 
This algorithm allows the designer to assign (mtr-1) eigenvalues of 
a system provided that it fulfills the sufficient conditions of the 
proposition presented in Chapter Three. Since the evaporator model 
fulfills these conditions, it was possible to assign all of the five 
of its eigenvalues. 

The various steps involved in the application of the method 
to the evaporator model have been described in Table 4.15. Vectors 
8) and B> refer to the column vectors involved in Step 1 and 
Step 2 of the algorithm and the asterisk superscripts designate 
those eigenvalues whose modes have been made uncontrollable in Step 2 
of the algorithm, (i.e., these eigenvalues were assigned in Step 1 
and did not change in Step 2). 

Initial attempts to apply the algorithm to the evaporator 
model encountered the difficulty of generating high gains in the 
controller matrices as illustrated by Runs 25 and 27. The controller 


matrix G shown in Table 4.18 contains gain elements which are 


ee) 
unacceptably large for the satisfactory control of the actual evaporator. 
As noted by Fallside and Seraji [13] the magnitude of the gain elements 
in the controller matrix depends on the ratio of the elements in the 
gulvectors IeSincelthe choice of g is arbitrary as long as the 
corresponding equivalent system remains controllable, this design 


freedom can be used to reduce the magnitude of the gain elements. The 


controller matrix obtained in Run 26 has resulted from an extensive 
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trial and error search for such a gi: Although most of the gain 
elements of Run 26 are considerably smaller than those of Run 25, one 
of them is still greater in absolute value than the limit of 100.00 
set in Section 4.2. 

Rather than assigning m of the sytem eigenvalues simulta- 
aneously in Step 1, they can be assigned recursively, one or more at 
a time. The question then arises as to how &) should be chosen for 
each recursive part of Step 1. In dyadic state feedback methods, where 
pairing of closed-loop and open-loop eigenvalues is possible, Simon 
and Mitter's [14] strategy of choosing g so as to maximize the con- 
trollability of the single eigenvalue to be assigned can be employed 
(cf., Section 2.4). However, in output feedback techniques one does 
not have a priori knowledge about which open-loop eigenvalue will be 
shifted to the desired closed-loop location. Thus, a modification of 
Simon and Mitter's approach [14] is needed. 

The unique strategy adopted in this thesis was to perform 
Step 1 of the algorithm developed in Chapter Three such that only one 
eigenvalue would be assigned recursively. Thus, rather than assigning 
three eigenvalues their desired closed-loop locations only one was 
assigned this value while the other two were assigned their respective 
open-loop locations. The three possible g's fulfilling equation 
(2.55) were then tested and the one giving the smallest gains was used 
in the final design. It should be noted though that in most of the 
runs only two recursive steps were performed in Step 1, since the 
resulting eigenvalue locations were deemed satisfactory. 

Run 27 and Run 28 clearly indicate the success of this 


approach in reducing the magnitude of the controller gains. [In fact, 
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in all of the runs where this approach was used, none of the controller 
matrix elements exceeded a value of 45 for the desired eigenvalues 
under consideration. However smaller desired eigenvalues resulted in 
larger gains, as would be expected. 

The major design options in this eigenvalue assignment method 
are: 

1) Choice of the desired closed-loop eigenvalue locations. 

2) Choice of the g's to be used. 

3) Choice of the number of recursive substeps to be used 
in Step 1 of the algorithm. 

The second design option has been utilized in order to 
ensure gain elements of reasonable magnitude in the controller matrices 
as described above. 

Runs 29-31 illustrated the fact that both the choice of 
desired closed-loop eigenvalues and the number of recursive steps 
used affect the closed-loop transient response characteristics con- 
siderably. This is in agreement with the arguments put forward in 
Chapter Two: the response characteristics of a system are not only 
governed by its eigenvalues but also by its eigenvectors whose direc- 
tions depend on the type of eigenvalue assignment scheme employed. 
It is obvious that the choice of the g's directly affects the closed- 
loop eigenvectors, but the analytical relation between these variables 
is not known, yet. 

Runs 32-35 are presented to indicate that even comparatively 
small changes in one of the five closed-loop eigenvalues of the system 
can result in considerably different response characteristics as shown 


in Figures 4.19 and 4.20. 
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Runs 36 and 37 depicted in Figures 4.21 - 4.24 are the best 
results obtained in a large number of applications (> 100) of the 
new algorith to the evaporator model. The reason for their inclusion 
in this thesis is to convince the reader that eigenvalue assignment 
can give excellent results if one is prepared to make exhaustive use 
of the available design options. 

Runs 38 and 39 use only the first part of Runs 36 and 37, 
respectively. Figures 4.25 and 4.26 clearly indicate the need of 
using a recursive scheme in this particular application where repeated 
eigenvalues with distinct eigenvectors occur. 

The following observations are made concerning the application 
of the eigenvalue assignment technique developed in Chapter Three to 
the evaporator model: 

1) Significant variations in the controller matrices and 
the closed-loop transient responses were observed. 

2) CF disturbances were more severe than TF and F dis- 
turbances and generally caused significant offsets in the controlled 
variables. This was due to the fact that in every run the steam flow- 
rate reached a steady state below its normal steady state value rather 
than above it. Bl and B2 were always properly manipulated. 

3) By changing the value of a small set of desired eigen- 
values, one could always reach a satisfactory closed-loop performance. 

4) Use of a greater number of linearly independent g's in 
recursive steps in general, tended to equalize the magnitude of gain 
elements. 

5) The large gain elements used in Run 29 caused some 


overshoot to the feedflow disturbance. 
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The simulation studies of this section clearly indicate that 
the objective of gaining some control over the entire set of open-loop 
eigenvalues is a more desirable one than the exact assignment of some 
of the eigenvalues. This point may be appreciated better in applications 
where the dimension of the system model is much larger than the number 
of available measurements and controls. In such instances the com- 
bination of exact and approximate eigenvalue assignment methods by the 
application of the following procedure may be useful: 

1) Assign max (m,r) of the system eigenvalues by the 
application of Step 1 of the algorithm described in Chapter Three. 

2) Determine B or k, by the application of Step 2 of 
the same algorithm in order to protect max (m-l, r-l) of the already 
assigned eigenvalues. 

3) Apply Jameson's [15] approximate eigenvalue assignment 
technique by employing Bo designed above. This Boar on cn will 
minimize the objective funetion defined in (2.73) and approximately 
assign those of the remaining max (n-mt+l, n-rtl) system eigenvalues 


which belong to both controllable and observable modes. 


4.8 Conclusions 

In this section some of the more important observations made 
in the simulation studies will be cited. Their interpretation and 
correlation will be presented in Chapter Five. 

Controllers derived via the application of Rosenbrock's 
approximate modal control method [8] produced closed-loop systems with 
similar response characteristics for any set of design options and 


disturbances considered. The location of the desired closed-loop 
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TABLE 4.6 


Design Specifications for Rosenbrock's Approximate 


Modal Control Method [8] 


Run 


Number Ordering of Eigenvectors* Desired Eigenvalues 
Wig e685! Hae ie 0.3 0.5 7 
a? CeCe eam Gener f 0.3 On 0.7 
3 Wa, Ee, hy igs 8 Oo5 0.7 
4 Wy Ms; M &, we 0.3 Ue OF 
5 Be Meee lee Ors a) Ore 
6 Ww, Ws Wo A) as Oi Oi 0.7 
: Wy My Wah Me Ue 0.00006 0.0003 0.001 
8 WW Wy We wy Q700006 0.00034. 0.001 
9 wg op a ae 0.8 ete) 0.8 
10 Wig aoa > Fy ale 0.5 5 Ons 
oe WW Wy MW, We 0.05 0.05 0.05 
12 Wie Np? We UH, MMe 0.00005 0.00005 0.00005 
13 W M W, WW, Oh 3! 0.4 0.6 
14 WW Wy Wy We Ort 0.4 On 
15 W YW Wy Ww WM, 0.4 0.0001 0.7 
16 W, Wy Wy Wy We 0.4 0.0001 0.0007 
17 Wy Moe hee 0.0004 0.0001 0.0007 
18 W W Wy Wy Ms 0.0001 0.4 O07 
19 W, MW Wy Ww, We 0.4 0.7 0.0001 
“The eigenvectors Rou emner Me eer the corresponding columns of 


matrix W in Table 4.4. 
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TABLE 4.7 


Response Characteristics for Rosenbrock's Approximate 


Modal Control Method [8] 


Run Corresponding Offset.in..C2' Offset in C2' 
Number Figures foro Disturbance, for CF Disturbance 
1 Ge a4 -0.013 0.014 
2 -0.018 0.020 
3 -0.018 0.020 
4 -0.013 0.014 
5 -0.030 0.033 
6 4.3, 4.4 -0.030 O}..033 
7 -0.009 0.010 
8 -0.009 0.010 
g AD, 4.6, Gat -0.043 0.048 
10 455, 426, 2578 -0.018 0.020 
8 and 4.9 
Ge Lge GalOue aed -0.009 0.010 
ud -0.009 0.010 
13 GG Cee -0.022 0.025 
14 -0.030 0.033 
15 WaOe OAc) -0.030 0.033 
16 -0.009 0.010 
17 -0.009 0.010 
18 -0.030 0.033 


19 BLOG a Le -0.009 0.010 
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Actual Closed-loop Eigenvalues for Rosenbrock's [8] 


Approximate Modal Control Method 
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Controllers for Rosenbrock's Approximate 
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Closed-loop Eigenvalues and Right Eigenvectors for Rosenbrock's 


Approximate Modal Control Method [8] 


Run 1 - Eigenvalues: 


0.502 


0.607 


TABLE 4.10 


0.65740.2691 


0.897 


Corresponding Right Eigenvectors (columnwise): 
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Table 4.10 (continued) 


Run 11 - Eigenvalues: 
0.050 0.050 0.62540.431i 0.897 


Corresponding Right Eigenvectors (columnwise): 
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TABLE 4.11 


Modal Control Method [9] 
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Design Specifications for Takahashi's Approximate 
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TABLE 4.12 


Response Characteristics for Takahashi's Approximate 


Modal Control Method [9] 


Run Corresponding Offset in C2' Offset in? C2” 
Number Figures for F Disturbance for CF Disturbance 

20 Gel, anes =O 2167 0.087 

20 -0.162 0.085 

22 G12.) Bais -0.091 0.049 
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24 aa Aaa -0.018 0.009 


ct) Too 
ry 
é 
U 


‘te bode 


‘SD. ar 388220 


tox: 


s5csG sua bid 7 ae 


Tal .O- 


Sdf.0- 


Te20, 0- 


SL .0= 
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Modal Control Method [9] 


Run Number 
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CHAPTER FIVE 


CONCLUSIONS 


5.1 Interpretation of the Simulation Results 
The different response characteristics observed in the 
application of eigenvalue assignment methods and Rosenbrock's 
approximate modal control method may be explained by the following 
arguments. 
Consider the dyadic representation of the open-loop system 
n 


matrix in terms of its eigenproperties, namely, 9 = ) hy i geN 5 
i=l 


The recursive application of the eigenvalue assignment method of 


Chapter Three to assign &% = min (m,r) system eigenvalues will 


result in the following closed-loop system matrix, 8 or, 4 
Heierkc PLAS CB ERNE HN Pe A Sob y « Gap 
CL {=1 ited i: t=. 2 i 


The directions of the {k,} in (5.1) depend on the set of desired 
closed-loop eigenvalues, whereas the ae can be assigned almost 
any set of elements. Thus, the closed-loop eigenvectors will, in 
general, be different from the open-loop ones and will highly depend 
on the choice of desired eigenvalues and the ig}, 

Application of Rosenbrock's approximate modal control 
method for the same purpose, on the other hand, will result in the 


following closed-loop system matrix (cf., (2.61)): 
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or 

n & ms ‘ 
bor = zi dy WP Sv, + ) k. h><f, C522) 
i=l] i=l 
: th ak -1 S 
where hy represents the i column of (vo By) sand f, repre- 
h — nan “a 
sents the i° row of (C Wy) hy hy and i, do not depend on the 


set of desired closed-loop eigenvalues, but Kk does. Therefore, 
for each set of desired closed-loop eigenvalues, in general, a dif- 
ferent set of closed-loop eigenvectors will be obtained. But, as 

a4 (l= 1525. ..,2) sapproacnes zero, ke Cia ohare Co). OWE 
attain the value of the open-loop eigenvalues which are being shifted 
(cf., 2.32). Thus, the second summation in (5.2), namely the con- 
troller matrix, will become a constant matrix which in turn means 
that the closed-loop system eigenvectors will assume fixed directions 
in the state space. Thus, in Rosenbrock's approximate modal control 
method [8], the sensitivity of the eigenvectors to the specification 
of the desired closed-loop eigenvalues diminishes when ai 
(i = 1,...,2) approaches zero and, thus uniformly changing response 
characteristics which approach an asymptotic limit are observed 

(cf., Runs 9-12). 

These observations cannot be made in the application of 
eigenvalue assignment techniques where the magnitude of the elements 
in {k, J of (5.1) may become unbounded as the desired eigenvalues 
are made smaller. Even for the same specification of the {g.} 
the eigenvector directions will be, in general, highly affected by the 
choice of the desired closed-loop eigenvalues, and thus significantly 


varying response characteristics can result, (cf., Runs 32-35). 


Similarly, since the eigenvector directions highly depend on the 


Bir i 


-Sroisten me: 
~~ EL ¢ t 
cae a ; a 
a i 1 r = 2% af viniie : * 
fiw (&...,8,f = 3) .S oven series tthe 
po - 
_ 


93341 sited sysp sisttw esutavnagis gool-=nago 3 


ae 1 ‘ ~ rrr. ¢ 
-roo ats vierspe . 4.2) ai solstnimpe tae 


i ey : ee z <4 ~ Vie 
sxud ol notdw xiatea tnesanes a st . exis ‘I H0 


BSULBV TE gf « 
' ; an 
a Se oe codanorage baa 
- ti ar 
io © ¢ a 4 ra ae. mG 4 _“! tS dase 1gga dotviw aak: 


ay, aan 
+ of sbam otf Jonnss seolva wisedo Saen we 
a a 


oe 


a3 > SbUI FaSsm srt 2848 esupianos3 eee iis — 

ulevasgis bsttesh aft es hshayodiy smoosd Yam (1. ye {3 # 

\,3) 3 To ookisorilssee sie ef4 tol neva ii aa S Df 
vides . janez of .sd Iitw enotstsex1tb Ag ov 


7 i” 
trspitiovnte ead & ,29ulsvoeests qeol-beeols bertiash efi2 ey 99. 
. . ae 7 
.(ct~St ane ,,i9) ,ifvess nao eslieltstjobreds sepegas } Be 
> 


, ee 
63 so basgeb Yiigia enbijzepth tososvasgis sia sante —ls 
To 
‘9 
. 


- 


129 


specification of the desired closed-loop eigenvalues, different response 
characteristics for CF and F disturbances are observed in different 
runs, i.e., for some runs, CF disturbances cause larger offsets and in 
other runs F disturbances are more severe. 

This then forms an explanation of some of the phenomena as 


could not be accounted for in Chapter Four. 


5.2 Comparison of Design Methods 


A rather detailed discussion and comparison of the various 
modal control and eigenvalue assignment methods considered in this 
thesis have been provided in Chapter Two. Thus, this section will only 
summarize the experience gained with the four methods simulated in 
Chapter Four. 

Both modal control and eigenvalue assignment methods increased 
the degree of stability and the speed of response of the system under 
consideration. The design options and design parameters involved in 
eigenvalue assignment methods were found to have a larger effect on 
the shape of the closed-loop system response than those involved in 
modal control methods. But currently there is no systematic way of 
exploiting these design options and parameters in order to improve 
the shape of response. 

Application of Rosenbrock's approximate modal control method 
resulted in favorable response characteristics provided that the 
desired real eigenvalues had reasonably small magnitudes. The amount 
of offset in the controiled variables decreased as the real parts of 
the closed-loop eigenvalues decreased. The closed-loop eigenvectors 


tended to assume fixed directions in the state space as the desired 
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real eigenvalues tended towards zero. Thus, the shape of response 
could not be affected to any great extent. The design options con- 
Sisting of the ordering of the eigenvectors in the eigenvector matrices 
and the pairing of closed-loop and open-loop eigenvalues influenced 
the response characteristics of the system, but to a smaller extent 
than the location of the closed-loop eigenvalues. The gains in the 
controller matrices had the physically expected signs and magnitudes, 
were always reasonably small, and did not vary significantly for the 
various options considered. Different disturbances had a similar 
effect on the closed-loop response characteristics of the system for 
the same control law. The successful application of Rosenbrock's 
approximate modal control technique in this case study was attributed 
to some of the desirable features of the open-loop eigenproperties. 
Application of the eigenvalue assignment algorithm of 
Chapter Three resulted in a wide spectrum of closed-loop response 
characteristics. The closed-loop eigenvector directions and the shape 
of responses were quite sensitive to the existing design options and 
design parameters. Thus, extensive utilization of these design options 
and parameters resulted in satisfactory controllers. The eigenvector 
directions and the response characteristics of the closed-loop system 
could not be correlated with the various design options and design 
parameters. The controller gains required for the assignment of a 
set of eigenvalues were always larger in absolute value than those 
resulting from the application of Rosenbrock's approximate modal con- 
trol method. Recursive application of the eigenvalue assignment methods 
was found to reduce the controller gains and to improve the response 


characteristics. Some suggestions were made in Chapter Four for the 
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successful application of this procedure. Different disturbances had 
a vastly different effect on the closed-loop response characteristics 
of the system. 

Application of Takahashi's approximate modal control method 
gave large offsets in the controlled variables which could not be 
reduced greatly by the specification of smaller real desired closed- 
loop eigenvalues. The combination of a mode-based model reduction 
technique with the ideal modal control method, on the other hand, gave 
excellent results. The generality of the results noted in this para- 
graph is rather doubtful, though, since extensive studies were not 
performed to evaluate these last two approaches. 

It is worth noting that approximate modal control techniques, 
in general, cannot guarantee the assignment of a set of eigenvalues. 
In fact,their application to some systems may even result in a less 
desirable set of closed-loop eigenvalues than the open-loop ones. 
Their success in this case study is due to the favorable eigenproperties 
of the open-loop system. By contrast the eigenvalue assignment method 
of Chapter Three guarantees the assignment of at least a subset of the 
system eigenvalues. Thus, its applicability is greater, and as demon- 
strated in Chapter Four it can provide excellent results if the avail- 
able design options and parameters are properly exploited and the 
unassigned closed-loop eigenvalues do not attain undesirable locations 
in the complex plane. 

The computational efforts involved in designing controllers 
is relatively small in both modal control and eigenvalue assignment 


methods. 
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5.3 Recommendations 

The simulation studies in this thesis involve a relatively 
low order model. Therefore, it has not been possible to demonstrate 
the usefulness of modal analysis in configurating the control system. 
Similarly, it has not been possible to fully demonstrate the benefits 
of the new eigenvalue assignment method which in most instances 
guarantees the assignment of a significantly larger number of eigen- 
values than in previously developed methods. Future work aiming to 
assess the real benefits of modal control and eigenvalue assignment 
should involve systems with higher order state space models. 

The greater importance of the closed-loop eigenvector direc- 
tions on the closed-loop response characteristics of the system be- 
came readily apparent during this investigation. It was realized that 
eigenvalue assignment did not provide satisfactory response character- 
istics unless extensive tuning was performed. Thus, systematic methods 
to affect the modes of the system rather than just its eigenvalues are 
greatly needed. The potential usefulness of the design freedom avail- 
able in choosing the g vector in fulfilling this objective is undeni- 
able. Preservation of the open-loop eigenvector directions in the 
closed-loop system is an arbitrary and unjustified objective unless 
the open-loop eigenvectors already have desirable directions. 

Exact eigenvalue assignment has been found to be an unnec- 
essarily ambitious goal. Rather than exactly assigning a few eigen- 
values, future work should attempt to gain some control over all of 
the system eigenvalues. 

Studies concerning determination of desirable closed-loop 


eigenvector directions for a specific system to fulfill such objectives 
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as disturbance rejection, high integrity, and low sensitivity to con- 


troller gain variations may prove to be very fruitful. 


NOMENCLATURE 


Alphabetic 

a coefficient in a characteristic polynomial 
A system matrix of a continuous-time model 

b coefficient in a characteristic polynomial 
b control coefficient vector 

B control coefficient matrix 

Bl first effect bottoms flow 

B2 second effect bottoms flow 

ic output coefficient vector 

C output coefficient matrix 

cl first effect concentration 

G2 second effect concentration 

CF feed concentration 

d disturbance vector 

D disturbance coefficient matrix 

S vector,defined,in (3.18) 

E matrix defined in (3.25) 

£ column vector of F 

F feed concentration 

5 modal output coefficient matrix 


F( ) matrix function defined in (3.5) 


ER dyadic matrix 
g column vector of a dyadic controller 
G controller matrix 
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h column vector of H 
H modal control coefficient matrix 
H1 first effect enthalpy 


HF feed enthalpy 


J, objective function defined in (2.72) 
J, objective function defined in (2.73) 
k controller gain 

k controller vector 

K modal domain controller matrix 

Q number of eigenvalues assigned, number of modes controlled 
L weighting on eigenvalues 

™m number of outputs 

m mode vector 

M fees on controller gains 

n number of states 

p number of disturbances 

Pp column of matrix P 

P analyzer matrix 


q( ) open-loop characteristic polynomial 


q vector defined in (2.26) 
Q controllability matrix 
r number of controls 


r( ) closed-loop characteristic polynomial 
ae column of R matrix 


R synthesizer matrix 
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modal control element 
modal control vector 


observability matrix 


modal disturbance element, time 


modal disturbance vector 
discretization interval 

control element 

control vector 

matrix defined in (3.17) 
left eigenvector 

left eigenvector matrix 

right eigenvector 

right eigenvector matrix 
State variable 

State vector 

output element 

output vector 

modal state variable 


modal state vector 
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Greek 

0, Scalar defined in (2.40), objective function defined in (2.70) 
B scalar defined in (2.51), objective function defined in (2.71) 
one Kronecker delta 

A partition of a matrix 

A discrete control coefficient matrix 

Q discrete disturbance coefficient matrix 

d discrete state coefficient matrix 

r eigenvalue 

A eigenvalue matrix 

Vv index 

E modal activation, scalar defined in (2.52) 

II product 

) sum 

t time 

Uy) degree of controllability 

Superscript 

: time derivative 

: perturbation variable 

* vector projection, arbitrary vector or nants 
T matrix transpose 

A first parcition of a matrix 

M second partition of a matrix 


-1 matrix inversion 
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Subscript 

_ vector 

= matrix 

c controllability 

d desired 

i element counter, run counter 
j element counter 

QL size of matrix partition 

fe) observability 

Symbols 


< = dor product 


>< dyadic product 
Abbreviations 

min. minimum 

max maximum 


sign signum 
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